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1.0 INTRODUCTION 

The purpose of this document is to present statistical characteristics of 
the SKYLAB S-193 altimeter altitude noise measurement process in a concise format for 
use and analysis by the Scientific Community. The data contained in this report 
is obtained from all three SKYLAB missions (SL-2, SL-3, SL-4). Section 2 
describes the statistical processing techniques, Section 3 presents summary 
residual statistics for 184 data arcs and Section 4 contains detailed noise 
statistics for 27 selected data arcs. Appendices A and B contain mathematical 
descriptions of the computer programs used for this data processing. 
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2.0 ALTIMETER DATA PROCESSING 


2.1 Measurement Description 

The altimeter range residuals are defined as the differences between the 
eight per second altimeter measurements to the sea surface and the calculated 
altitude from the spacecraft to a reference ellipsoid. This calculated altitude 
is derived from a short arc orbit determination xislng S-Band and C-Band tracking 
data as described in reference [2]. The calculated altitude Includes corrections 
for pulsewidth, bandwidth and other operating parameters. A 2.79 meter correction 
for refraction was also applied. Errors in the calculated altitude are constant 
over the altimeter data arc thus contributing a constant bias error. The noise 
on the altimeter measurements is a function of received signal to noise ratio, 
the statistical Independence of the returns and tracking servo bandwidth. For a 
near beamwidth limited system such as the S-193 (in the 100 ns pulsewldth mode) 
the received peak power is a very sensitive function of pointing angle with 
respect to nadir. An increase in pointing angle will (1) introduce a bias error 
(2) reduce signal power and (3) change the return signal waveshape. The reduction 
in signal power will increase the measurement noise variance and reduce the track- 
ing servo bandwidth thus altering both the noise magnitude and power spectrum. 
However, the change in the received waveshape will have the largest effect on the 
servo bandwidth for it changes the characteristics of the error detection. 

Altimeter range measurements were analyzed for three major operating modes, 
(Table 2.1) modes 1, 3 and 3. Each operating mode has three data acquisition 
submodes (DAS-1, DAS-2 and DAS-3) which vary in length from 16 seconds to 103 
seconds. In this report a data arc is defined as the residuals recorded during 
one altimeter mode/submode. Each such arc is uniquely identified by a map number 
as assigned in reference [2] and the submode designation. Hardware problems which 
influenced the analysis and statistical summaries presented in this report in- 
cluded the following: 

• The pulse compression (Mode-5 DAS-2) was inoperative for all data 

arcs prior to Map 73. As a result for these data arcs the altimeter 
processor operating characteristics were nearly identical to Mode-1 
DAS-2 

The antenna gain was reduced to 28 db during SKYLAB IV. This affected 
every data arc after Map 76. 

A more detailed description of the altimeter experiment and hardware is given 
in reference [2J. 
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In Mode 1 DAS-3 the antenna pointing angle is offset .431® in pitch from its 
nadir pointing position to evaluate the effect of known pointing angle variation. 
Pointing angle can be estimated (1) from an analysis of the tracker response and 
automatic gain control (AGC) and (2) from an analysis of the return waveform. The 
pointing angles as tabulated in this report, have been estimated by the first 
method which is accurate to within a few tenths of a degree for pointing angles 
from .75® through 1.7®. 


2.2 Evaluation of Residual Statistics 


In addition to broadband (Bandwidth 2 Hz) measurement noise the altliseter 
residual data contains low frequency effects due to geold variations. Thus 
evaluation of the noise statistics consisted of (1) selection of a reasonable 
nund>er of data arcs (approximately 184) • which visually contain a minimal 
amount of geold variations, (2) computation of stmunary residual statistics for 
all of these data arcs, and (3) hl^ pass filtering and computation of detailed 
noise characteristics for twenty-seven typical data arcs spanning all nine altim- 
eter operating modes/submodes and three pointing angle ranges (low, moderate, 
and hl^). 

2.2.1 Sunanary Residual Statistics 

The sumnary residual statistics as compiled from the RAP program output 
include: 

* Mean 

• Standard deviation 

* F Statistic for assessment of variance homogeneity 

• Chi square statistic to assess distribution normality 

• Listing of significant power spectral.«omponents 

* Listing of significant autocorrelation coefficients 

A detailed description of the computation as performed In RAP Is contained In 
Appendix A. The F statistic is computed by partitioning the data arc into four 
segments, computing a mean and sum of squares about the mean for each segment and 
then conq>aring the pooled sums of squares about the segment means with the sum 
of the squared deviations of the segment means from the overall means. The chi square 
statistic measures the goodness of fit of the residual histogram to the normal 
distribution. 

In RAP the power spectrum Is computed and scaled such that if the residuals 
are a white noise process then each spectral component Is an independently dis- 
tributed chi square variate with two degrees of freedom. The raw power spectrum 
(l.e. before application of any window functions) is then divided into fifty cells 
containing an equal number of spectral coefficients. The coefficients are susned 
within each cell. Under the white noise assumption these sums are also Indepen- 
dently distributed chi square variates with degrees of freedom equal to the sum of the 
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degrees of freedom of the coefficients within that cell. Thus if p coefficients 
are summed within a cell the resultant sum would have 2p degrees of freedom. 
Significant frequency bands are then defined as those cells where this sum ex- 
ceeds the 99th percentile of the appropriate chi square distribution. Thus, if 
the observed resldiial time series is not a white noise process, this approach 
identifies those frequency bands where the signal power is concentrated so as 
to cause significant deviations from the white noise assumption. 

Similarly, for a white noise process the autocorrelation coefficients are 
asymtotically normal with mean zero and standard deviation l/vfT'where N is the 
sample size. Significant autocorrelation coefficients are defined as those 
coefficients whose absolute value exceeds 2 /^ (i.e. two sigma). As with the 
spectrum, this approach simply identifies the coefficients which are significantly 
different from their expected value (zero) under the white noise assumption. 


These summary statistics are presented and discussed in Section 3. 



2.2.2 Detail Noise Statistics 


To evaluate S-193 measurement noise characteristics In addition to computa- 
tion of the above statistics, the longer data arcs (i.e. longer than 24 seconds) 
subject to detail analysis were high pass filtered bv the program SAMPFL (Appendix 
B) to remove any low frequency geold effects. This high pass filtering Is 
accomplished by SAMPFL by (1) low pass filtering the data with a 99 point 
quadratic midpoint filter and (2) computing the high pass ovtput as the difference 
between the input data and the low pass output. The low pass power transfer 
function plot (Figure 2.1) shows that cut off frequency for this filter is .0845 Hz. 
In reference [2] a low pass filter with the half power frequency of .164 Hz 
is recomnended for recovery of geold features from the SKYLAB S-193 data. Since 
the purpose of this study was to evaluate the S-193 measurement noise statistics, 
the data arcs processed had no prominent geold features, therefore the .0845 Hz 
cutoff was adequate for rejection of any small geold features or trends and yet 
retained most of the measurement noise power. 


In reference (2] the S-193 tracking loop frequency domain transfer function 
is given as 


H,,(jw) 


701W+280 sin v^/2 
8(jw)^+81jwf280 "TTT” 


^-jwT/2 


( 2 . 1 ) 


where w is the radian frequency (2irf) and T is the sampling period of .128 
seconds. From the plot of |H^(Jw)| (Figure 2.2) it is obvious that the 99 point 
highpass filter removes only a small fraction of the broadband noise power. 


After high pass filtering (where applicable) the following were computed 
and displayed. 

* Summary statistics as defined in Section 2.2.1 

* Time series plot 

« Power spectrum plot 
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Autocorrelation function plot 


• Autoregressive model power transfer function plot 

The summary statistics included mean, standard deviationp F statistic, chi 
square statistic, significant power spectral coefficients and significant auto- 
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POWER TRANSFER FUNCTION OF 99 POINT QUADRATIC MIDPOINT FILTER 

FIGURE 2. 1 
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FIGURE 2.2 

POWER TRANSFER FUNCTION 
S-193 PROCESSOR 
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correlation coefficients as defined above in Section 2.2.1. C/.LCOMP plots of 
the high pass filtered (where applicable) time series, power spectrum and auto- 
correlation function were generated using the plot tape output option of the 
RAP program. For plotting purposes, to present a smoother graph, the raw power 
spectrum was modif J d by two applications of the Hanning window to the fourier 
coefficients as described in Appendix A, Section 2.6. 

To provide an analytical estimate of the noise power spectra, an 
autoregressive model was fitted to the residual time series. 

The power spectral density requires 1 + N/2 coefficients (where N is the sample 
size) to describe the detailed frequency characteristics of the data and the 
autocorrelation function requires N-1 coefficients to describe the time domain 
correlation structure. These are descriptive statistics which make no assump- 
tion about the underlying process which generates the time series. The auto- 
regressive model [ 1 ] assumes that a stationary time series {X^} can be represented 
as finite sum of the form 


'^t ■ ^l*t-i 


^2*t-2 


p t-p 


( 2 . 2 ) 


where X in the observed value at time t, c is a zero mean white noise process 
^2 *■ 

with variance a and the { 41 ^; l*l, 2 **p} are called the autoregressive coefficients. 
This is equivalent to assuming that {X^} is generated by passing a discrete 
white noise process {e^} through a digital filter with the transfer function 
H(Z) given by 

H(Z) — ir — (2.3) 

• • •— (J) z 

1 2 p 

where Z ■ , f is the frequency in Hz, T is the sample time in seconds, and 

j - /-I. 


A linear least squares procedure is implemented in RAP (Appendix A) to 

estimate the autoregression coefficients from the computed autocorrelation 

coefficients. For purposes of this study it was arbitrarily decided to estimate 

the first twenty autoregression coefficients. From these estimated coefficients 

the power transfer function |h(Z)| was computed and plotted over the frequency 

range (0-3.9 Hz). This plotted power transfer function is related to the power 

2 

spectrum of the digital filter output, N(f), by the scale factor 0 as 
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N(f) - 

These detailed noise statistics are 


2 - 



Box and Jenkins, Time Series Analysis. Forecasting and Control , Holden* Day, 
1970. 

McGoogan, Leltao and Hells, "Sunmary of SKYLAB S-193 Altimeter Results", 
NASA TM-X-69355, NASA Wallops Flight Center, February 1975. 


3.0 STATISTICAL SUMMARY 


This section presents statistical summaries of altimeter range residual 
data from 184 data arcs. Table 3.1 gives the number of analyzed arcs for each 
mode/submode. Within a given mode/submode the antenna off nadir pointing angle 
is the principal factor which determines the altimeter hardware signal processing 
parameters. Therefore, in this section the residual statistics are tabulated 
as a function of pointing angle with the pointing angle ranges defined as 
follows : 

Low • pointing angle $ .6 degrees 

Moderate - .6 degrees < pointing angle $ .9 degrees 

High - .9 degrees < pointing angle 
In all cases the tracker determined pointing angle is used. 

Tables 3.2 through 3.10 give the date (year, month, day), start time 
(hours, minutes, seconds), stop time, pointing 

angle and summary statistics as extracted from the RAP program output. The F 
statistic, which assesses the variance homogeneity over the data arc, has (3,N-3) 
degrees of freedom where N is the sample size. Large values of this F statistic 
would indicate a lack of such homogeneity. For the sample sizes given in tables 
3.2 through 3.10 the 99th percentile of the F distribution is approximately 3.78. 

F values greater than this are marked with an asterisk (*). Residual normality 
is assessed by computing a chi square statistic which measures the difference 
between the residual histogram and the normal distribution. The table entry 
(Normality Test) is the value of the chi square cumulative distribution function 
associated with the computed statistic. As discussed in Section 2.2.1, the 
significant frequency bands are those frequency bands where the grouped power 
spectral coefficient exceeds the 99th percentile of the chi square distribution 
with appropriate degrees of freedom. Significant autocorrelation coefficients 
are defined as those which exceed the plus or minus two sigma value for a white- 
noise process. If an autocorrelation function has all large positive coefficients 
the process is obviously non-statlonary. The table entry contains either the lag 
indices of the significant coefficients or the abbreviation NON-STA. indicating 
that the process is non-stationary. Usually this non-stationarity occurs when 
the data contains geoid features which are large in comparison to the measurement 
noise level. 



Section 3.1 presents the distribution of residual signas within each 
pointing angle range (low* moderate and high), and Section 3.2 describes the 
distribution of significant frequency bards. 


TABLE 3.1 

NUMBER OF DATA ARCS FOR 
EACH MOOE/SUBMOOE 


MODE 

POINTING 

ANGLE 

RANGE 

SUBMOOE 

DAS-1 

DAS-2 

DAS-3 

MODE 1 

LOW 

6 

6 

6 

MODERATE 

12 

10 

10 

HIGH 

5 
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2 

MODE 3 
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LOW 
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HIGH 
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MODE 5 

LOW 

10 
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TA8LF 3^ 

MODE 1 LOM POINTING ANGLE 



V 





NM> 

PASS 

MODE 

DATE 

START 

STOP 

POII^ 

(DECS) 

SAMP. 

SIZE 

^CAN 

(KTRS) 

SI6VA 

(MTRS) 

F 

STATISTIC 

NORM. 

TEST 

SIG. 

FREQ. 

BANOS 

(HZ) 

SIGNIFICANT 

AOTOCCRRELATION 

COEFFICIENTS 

, j 

i 

. 

1 

B 

1 D4S-1 

730604 

171U1 

171155 

D 

200 

46.7 

1.89 

2.33 

.979 

0-.0837 

fON ST A. 

1 

1 0AS*3 


171302 

171401 

■ 

400 

57.0 

1.73 

.731 

.361 

0-.059 
.391-. 449 

NON STA. 

36 

22 

1 OAS-I 

730902 

142954 

143034 

<.6* 

wSm 

9.17 

1.37 

1.16 

.641 

0-.0376 

10.11 


1 OAS-2 


143040 

143139 

352 

7.19 

1.21 

2.96 

.666 

fOONE 

8.19 


X D4S-3 


143157 

143243 


304 

8.76 

1.41 

.277 

.680 

0-.128 
.771-. 822 

NON STA. 

39 

22 

1 OAS-1 

730902 

143342 

143416 

<•6* 

208 

17.9 

1.95 

1.22 

.795 

0-.112 

NON STA. 

1 OAS-2 


143430 

143530 

408 

22.5 

1.20 

2.53 

.827 

0-.134 

NON STA. 

1 OAS-3 


143533 

143633 

408 

23.3 

1.30 

.344 

.746 

0-.0574 

NON STA. 

i i 

40 

22 

OAS-3 


143913 

144013 

<.S* 

408 

12.4 

1.36 

1.82 

.727 

0-.057 

NON STA. 

41 

22 

1 OAS-2 

730902 

144149 

144251 


416 

3.48 

2.02 

1.95 

.938 

0-.131 

NON STA. 

OAS-3 

144253 

144353 

<.5* 

408 

-5.1 

3.03 

.214 

.991 

0-.134 

NON STA. 

— ^ 

65 

55 

1 OAS-I 

731201 

173026 

173121 

<*5^ 

364 

31.0 

1.77 

.603 

.692 

0-.0429 

NON STA. 

DAS-2 

173127 

173210 


260 

30.5 

1.72 

2.07 

.526 

0-.0279 
.837-. 8651 
1.172-1.200 

1,339-1.367 

1.8.9.12. 

17,20 

66 

55 

i(2)0AS-t 

731201 

173356 

173442 


268 

12.7 

2.90 

.992 

.991 

0-.136 

NON STA. 

t<2)0AS-2 


173445 

173546 

<. 5 * 

416 

-3.16 

2.13' 

5.11* 

.956 

0-.056 
.300-. 357 
.826-. 883 
1,05-1.11 

NON STA. 


1 

< 

I 













































KftP 

PASS 


DATE 

START 

STOP 

POIOT 

(DEGS) 

86 

55 



173549 

173649 


91 

57 

1 OAS-1 

731202 

182231 

182256 


1 0A5-2 

182259 

182323 

<.5* 


TABLE 3.2 (CONT. ) 


SA.VP. 

SIZE 

4oa 


t92 


164 


MEAN 

(f>rrRS) 


-.03 


36.9 


37.9 


SIGMA 

itms) 


1.95 


1.84 


1.66 


statistic 


.549 


.226 


1.21 


NORM, 

TCST 


.945 


.989 


.7X8 


SIG. 

FREQ. 

BANOS 

(HZ> 


0*.057 
.536-. 594 
.609-. 747] 
.919-. 977 


.163-. 203 
.448 


• 425 
.934 


SIGNIFICANT 

AUrXCRRELATiaN 

COEFFICIENTS 


X.2.8-18 


1-5 


1.7.8.12 



1 
















TABLE 3.3 


MODE 1 MODERATE POINTING ANGLE 



PASS 

KODE 

DATE 

START 

STOP 



(NTRS) 

SIG^A 

r 

STATISTIC 


SIG. 

FREQ. 

B.VCS 

(HZ) 

SIGNlFICAVr 

AUTOCCRRELATIO': 

COEFFICIENTS 

201 

97 

1(2)0AS-1 

740131 

145526 

145610 

.75* 

240 

3.66 

2.12 

2.02 

.997 

0-.033 
.846-. 879 




X(2)0A3-2 


145614 

145714 


364 

B 

2.39 

B 

.575 

0-.1073 

NCN STA. 



K2)0AS-3 


145718 

14SS16 


342 

-12.2 

2.59 

3.11 

.303 

0-. 114 

NON STA. 

202 

97 

X(3)0A< -1 

740131 

150027 

150110 


260 

-25.7 

1.78 

.774 

.927 

0-.0279 

^ 1-4, 6-8. 

12-13.15- 
17.19 



K310AS-2 


150115 

! 

150215 

.75* 

360 

-27.6 

1.66 

! 

2.64 

.539 

0-.043 

NON STA. 



l(3)0AS-3 

1 

150216 

150318 


364 

-30 

1.66 

1.27 

.625 

! 

.193-. 236 
.645-. 687 

1,2,8.9.11. 

12 

203 

97 

l(A)OAS-l 

730131 

150921 

150505 


240 

-31.6 

1.66 

.243 

.997 

0-.0325 

1.107-1.139 

1.8.14 



l(A)OAS-2 


150509 

150609 

.75* 

408 

-35.5 

2.34 

1.83 

.766 

0-.057 

NON STA. 



K410AS-3 

. 


150611 

150711 


364 

-39.3 

1.62 

.961 

.961 

.258-. 300 
.451-. 494 
.708-. 816 
1.03-1,07 

■ 






























































































TABLE 3, 3 (CONT. ) 


MAP PASS 


l(S>OAS-2 


l(5)DAS-3 


l(6)OAS-2 


t(6)OAS-3 


206 9T 1(T)0AS-2 


l(7)0AS-3 


START 

STOP 

150808 

150851 

150856 

150955 

150959 

151058 

151157 

151249 

151252 

151345 

151347 

151448 

151630 

151725 

151734 

151834 


POIMT SAMP. MEAN 

(DECS) SIZE (MTRS) I (NTTRS) I STATISTIC 


2*0 -40.8 1.60 


.5 2.38 


SIG. 

FREQ. SIGNIFICANT 

BANOS AUTXORRELATION 
<H2) COEFFICIENTS 


352 -48.3 1.81 


280 -48.0 
.75* 308 -46.1 1,73 

360 -45.4 1.92 


.75* 312 -35.3 2.27 


.208 2.41-2.44 


.929 0-.137 

.859-. 918 


.880 O-.lll 

.400-. 443 
.533-. 644 
1.20-1.24 


.918 0-.106 

.382 -.488 


.581 0-.050 

.977-1.03 


.901 0-.043 

.586-. 629 


NON STA. 


1*2»12-14 


NONE 


1.2.17 


NON STA. 


NON STA. 


NON STA. 





























































































i 


f 


i 


i 


I 

00 


KAP 

PASS 

WOOE 

DATE 

START 

STOP 

207 

97 

l(3)DA5-l 

740131 

152200 

152252 

208 


1(9)0AS-1 

740131 

152553 

152637 



H9)D4S-2 


152640 

152741 



l(9)DAS-3 


152744 

. 

152343 

212 


1(13)DAS-1 

7430131 

154059 

154144 

213 

97 

l( I4}DAS-l| 

740131 

154443 

154527 



l(14)0AS-2 


154530 

1 

154630 



l(14)0AS-3 


154634 

154733 

214 


1(15)0AS-1 

740131 

154832 

154916 








l(lS>0AS-2 


154919 


issoia 


TABLE 3.3 (CXDNT. ) 


POINT 

(DECS) 

SA^5». 

SIZE 

^’^AN 

<:<TRS) 

SIG*<A 

(KTaS) 

F 

STATISTIC 

NORy. 

TEST 

SIG. 

FREQ. 

BATCS 

<H2) 

SIGNIFICANT 

AlTCCCFSELATICN 

CCEFFICIENTS 

.75* 

240 

21.72 

2.40 

.577 

.370 

0-.098 

NON STA. 

.75* 

240 

45.96 

1.67 

1.19 

.995 

NCNH 

UlO 

416 

43.1 

2.42 

2.26 

.992 

0-.056 

NCN ST A. 

400 

31.6 

3.48 

1.08 

.875 

0-.137 

NON 5TA. 

.75* 

294 

-60.1 

1.78. 

1.70 

.401 

0-.0266 
.79 7-. 824 

NON STA, 

i 

.75* 

240 

-33.0 

2.57 

1.04 

.651 

0-.162 

NON STA. 

360 

-23.1 

2.48 

.173 

.195 

0-.174 

1 

i 

NON STA. 

352 

-18.2 

2.35 

1.70 

.622 

1 

O-.lll 
.932-. 977 1 

1 

NON STA. 

• 

• 

cn 

• 

240 

2.03 

2.27 

.774 

.931 

.326-. 350 1 

.911-. 944 

1.14.18.19 

294 

1.30 

3.56 

2.62 

.624 

0-.132 
.21 3-. 292 
.425-. 452 
.531- .558 
.691-. 717 

NON STA. 











1 

SIGN!A 

(Naas) 

F j 

STATISTIC ! 

NCHM. 

TEST 

SIG. 

FREQ, 

BANOS 

<H2) 

SIG.NlFICArCT 

ALrrOCORRELATION 

CCEFFICIENTS 

■ 

4.60 * 


0- . 057 
. 153-.287 
.383-.4<.0 
.536-. 670 
.766-. S23 

I- 3, 5-7, 

II- 16 

2.53 

3.64 

.974 

.698-. 732 
1.12-1.50 

1- 

1,13,14 

2.2 

.916 

! 

1 

1 

.930 

^ 0-. 0434 
.391-. 434 
.521-. 629 
.701-. 825 

NON STA. 

2. IB 

1.02 

.999 

0-.114 
.342-. 386 
.430-. 731 

1,2,3,11-15 


. 022-. 868 















































TABLE 3.4 


MODE- 1 HIGH POINTING A.^ES 


2S 


33 


100 


PASS 


17 


19 


02 


MODE 


1 DAS-1 


1 OAS-2 


1 OAS-3 


1 OAS-2 


1 DAS-3 


I OAS-l 


b C»AS*2 


n DAS-1 


DATE 


730610 


7^0o<j9 


730S12 


731207 


START 


162619 


162901 


163006 


136806 


U 36656 


P23S02 


P23703 


166610 


STOP 


162857 


163002 


163106 


136667 


136966 


023567 


023606 


166653 


POINT 

(DEGS) 


1.2* 


1.3^ 


l.C“* 


1.3* 


Sam*, 

SIZE* 


192 


366 


352 


312 


600 


262 


666 


266 


^CAN 

(MTRS) 

69.0 


96.5 


166.3 


6.75 


57.3 


-69.6 


-30.9 


65.0 


SIGMA 

<KTRS) 

5.30 


6.23 


6.55 


5.60 


6.52 


3.03 


5.91 


7.69 


m 


i 


F 

STATISTIC 

NORM. 

TEST 

SIG. 

FREO. 

BANDS 

(HZ) 

SIGNIFICANT 

AUTOCORRELATION 

COEFFICIENTS 

1.16 

.796 

0-.081 
.607*. 529 
.651-. 692 
.695 

1-3. 6-6. 
10-13 

1.62 

.101 

0-.666 

NON STA. 

»0.13* 

.995 

0-.377 

ALL 

7.79* 

.999 

0-.500 

1-7 

12-17,19 

5.96* 

.207 

.0781-. 371 
.669-. 605 

1-5, 

.607 

.976 

.322-. 355 
.666-. 678 
.839-. 872 

ltl5 

5. 26* 

.761 

0-.45J 
.534-. S48 

1-6,11-19 






J 


3.95* 


981 


209-.653 
626-. 663 


1 - 6 , 6-12 































































3-12 


mm* 







> 

1 


■t 


TABLE 3.5 

MODE 3 LOW POINTING ANGLE 


r-w 

PASS 

KGDE 

DATE 

start 

STOP 

'^INT 

COEGS) 

SA?>P. 

SIZE 

MSAN 

(NTRS) 

SICXA 

{^^^R5) 

F 

statistic 

NORM. 

TEST 

SIG. 
FREQ. 
SANDS • 
(HZ) 

SIGNIFICANT 

AUTCCCRRELATIGN 

CCi£FFICIE\TS 

45 

24 

3{3)DAS-1 

730903 

154440 

154454 


104 

-81.6 

1.45 

.952 

.990 

NONE 

18 



3(3)045-2 


154455 

154517 

<.6* 

114 

-70.4 

1.57 

1.51 

. 161 

• 0685 

1,2. 4-8 
10.11,14. W 



3(3)CAS-3 


154600 

154703 


435 

1 

-es.7 

1.79 

4. 17* 

.991 

0-.054 

MON ST A. 

89 

57 

3 OAS-l 

731202 j 

161710 

161724 

<•5* 

1 

90 j 

[ 34.9 

1.71 

1.64 

.586 

.347 

14 



3 04S-2 


181726 

181748 


114 

1 48.6 

1.94 

2.91 

.939 

N&NE 

1 


TABLE 3.6 


MOOE-3 MOCCRATE POINTING ANGLE 



PASS 

MODE 

DATE 

START 

STOP 

POINT 

(DECS) 

SAW>. 

SIZE 

fCAN 

(MTRS) 

SIGMA 

(MTRS) 

F 

STATISTIC 

NORM. 

TEST 

SIG. 

FREQ. 

BANOS 

(HZ) 

SIGNIFICAMi 

AUTOCORRELATION 

COEFFICIENTS 

10 

e 

3 DAS-1 

730611 

152510,5 

152525 

.75* 

112 

“29,8 

1.44 

.581 

.704 

.349 

NONE 

3 0A5-2 

152526.5 

152549 

170 

“30.2 

1.42 

.621 

.838 

1.65 

1*14 

3 OAS-3 

152552.5 

152631 

198 

“31.7 

1.55 

.331 

.610 

0-.395 

6.8.10*11 

18.19 

m 

24 

3<l)0AS-3 

730903 

153623 

L53601 

.8* 

104 

“15.8 

4.45 

1,37 

1. 

0-.067 

NON STA. 

1 44 
— 

24 

3(2)0AS->2 

730903 

154022 

154043 

.8* 

no 

“39.5 

1.64 

1.46 

.506 

.071 

NON STA. 





3-14 


TABLE 3.7 


MOOE-3 HIGH POINTING ANGLE 


MAP 

PASS 

MODE 

DATE 

START 

STOP 

POIWT 

(DECS) 

SAMP, 

SIZE 

MEAN 

(MTRS) 

SIGMA 

(MTRS» 

F 

STATISTIC 

NORM. 

TEST 

SIG. 

FREQ. 

BANOS 

(HZ) 

SIGNIFICANT 

AUTOCORRELATION 

COEFFICIENTS 

8 

7 

3 OAS-1 

730610 

143629 

143639 


76 

130.8 

7.39 

4. 7* 

.997 

.206 

.411 

I. 2. 3. 4. 5, 

II. 12,13,14,15 

3 OAS-2 

143640 

143704.5 

1.45* 

184 

131 

9.2 

3.69 

.999 

.0849 
.170-. 382 

ALL 

1 

12 


3(1)DAS-1 

730612 

131001 

131015 


108 

50.4 

2.79 

.428 

.916 

.940 

1 



3(l)OAS-2 

131016 

131039 

1.1* 

176 

50.9 

2.25 

.100 

.451 

.133. .488 

1 


■ 

J(l)DAS-3 

131042 

131203 


526 

52.0 

2.69 

1.30 

.972 

0-.061 
.458-. 519 
.687-. 742 

1.11.12 

13 

1 

3(2)OAS-l 

730612 

131536.5 

131550 


104 

40.5 

4.27 

.541 

.845 

.751,. 977 

1.2.5 

3(2'o/- 2 


131551 

131614 

l.l* 

176 

40.7 

3.99 

1.57 

.300 

. 133,. 266, 
.400, .444. 
.666 

1.2.3 

3(2)DAS-3 

131621 

131V55 

612 

43.1 

6.45 

10,39* 

1.0 

0-.600 

1-6,14-20 

29 

17 

3 DAS-1 

7308CQ 

131512 

135125 

m 

100 

42.0 

3.00 

3.70 

.990 

.469 

1,17,18 

3 DAS-2 


131527 

135150 

176 

55.1 

2.91 

.402 

.986 

.355, .710 

1.6,7 

75 

40 

3 OAS-1 

730914 

170517 

170531 

* .0* 

104 

47.7 

5.40 

5.46* 

.984 

.0751,. 376 

1-5 

3 DAS-2 


170532 

170555 

176 

59.5 

6.83 

3.12 

.786 

.044-. 088 
.355-. 533 

1-4,6-10, 

14-20 

3 OAS-3 


1 70557 

1 70655 

256 

63,9 

6.06 

3.95* 

.960 

. 183-.458 

ALL 

135 

81 

3 CmS-2 
J 

740111 

172926 

172947 

1.0* 

160 

52.3 

2.35 

.331 

.893 

HI 

1.7,10 
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TABLE 3.7 


MOOe-3 HIGH POINTING ANGLE 


MAP 

PASS 

MODE 

DATE 

START 

STOP 

POIf^ 

(DECS) 

SAMP. 

SIZE 

MEAN 

(MTRS) 

SIGMA 

(MTRS) 

F 

STATISTIC 

NORM. 

TEST 

SI6. 

FREQ. 

BANDS 

(HZ) 

SIGNIFICANT 

AirrOCORRELATION 

COEFFICIENTS 

8 

7 

3 OAS-1 

730610 

143629 

143639 


76 

130.8 

7.39 

4.7* 

.997 

• 206 
.411 

I, 2. 3. 4. 5, 

II, 12,13,14,15 



3 OAS-2 


143640 

143704.5 

1.45* 

18* 

131 

9.2 

3.69 

.999 

.0849 
.170-. 382 

ALL 

12 

9 

3(1)0AS-1 

730612 

131001 

131015 


108 

50.4 



2.79 

.428 

.916 

.940 

1 



3( 1 )0AS-2 


131016 

131039 

l.l* 

176 

50.9 

2.25 

.100 

.451 

.133».488 

1 



3(l)DAS-3 


131042 

131203 


526 

52.0 

2.69 

1.30 

.972 

0-.061 
.458-. 519 
.687-. 742 

1811.12 

. 

13 


3(2)OAS->l 



131550 


104 

40.5 

4.27 

.541 

.845 

.751,. 977 

1.2.5 



3(2V^^^ 2 


131551 

131614 

l.l* 

176 1 

40.7 

3.99 ! 

1.57 

.300 

. 133, .266, 
.400, .444, 
.666 

1.2,3 



3(2)0AS-3 


131621 

131/55 


612 

*3.1 j 

6.45 

10.39* 

IB 

0-.600 

1-6,14-20 

29 

17 

3 OAS-1 



135125 

■ 

100 

42.0 

3.00 

3,70 

.990 

.469 

1,17.18 



3 OAS-2 



135150 


176 

55.1 

2.91 

.402 

.986 

.355, .710 

1.6,7 

75 

40 

3 OAS-1 

730914 

170517 

170531 


104 

47.7 

5.40 

5. *6* 

.984 

.0751,. 376 

1-5 



3 DAS-2 


170532 

170555 


176 

59.5 

6.83 

3.12 

.786 

.044-. 088 
.355-. 533 

1-4,6-10, 

14-20 



3 OAS-3 


170557 

1 70655 1 


256 

63,9 

6.06 

3.95* 

.960 

.183-. 458 

ALL 

135 

81 

3 Df\S-2 
J 

740111 



■ 

160 

52.3 

2.35 

.331 

.893 

■ 

1,7,10 



















































































TABLE 3.8 





r 

H* 


i 


MODE 5 LOW POINTING ANGLE 



PASS 

MOOc 

OATH 

START 

STOP 

POINT 

(OESS) 

SA^^». 

SIZE 

P£AN 

SIGMA 

<MTRS) 

F 

STATISTIC 

NORM. 

TEST 

SIG. 

FREG. 

BANOS 

(HZ) 

SIGNTFICAVr 

AUTOCORRELATION 

COEFFICIENTS 

It 

■ 

5 OAS-2 

730612 

130226.32 

130319 

■ 

352 

40.] 

1.30 

.226 

.993 

o-,in 

NON. STA. 


1 

5 OAS-3 


130336 

130426 

■ 

280 



14,6^ 

.908 

.in-. 140 

.335-. 363 
.614-. 698 

1.2.5.6.10 

33 

17 

5 OAS-1 

730809 

135828 

135843 

■ 

114 

17.3 

2.37 

.995 

.497 

.548 

1.13 



•r DAS-2 


135845 

140025 

B 

672 

6.07 

2.99 

til 

• 

.999 

0-.128 

NON STA. 

32 

18 

5 OAS-3 

73081 1 


153941 

B 

184 

-37.5 

.810 

5.56* 

1 

.425 

.977 

1.10 

1.4.11, 

16,17 

35 

21 

5(1)DAS-1 

730901 

152247 

152302 

ra 

104 

16.7 

1.22 

1.27 



N0N2 

46 

25 

5(l)OAS-l 

730904 

145138 

145152 

D 

108 ; 

-7.86 

1. 11 

1.16 

.430 

.0723-. 145 

1.2. 5. 6, 7 



5(UDAS-3 


15S337 

155428 

i 

■ 

342 

-29.1 

1.65 

2.60 

1.0 

0-.114 

NON STA. 

«7 

25 

S(2)OAS-l 

730904 

1455U 

145523 

■ 

i 

92 j 

-32.3 

I.IO 

.768 



16 



S(2)DAS-2 


145526 

145706 

B 

704 

-39.6 

3.04 

6. 98* 

1.0 


NON STA. 


■ 

5(2)DAS-3 


145709 

145759 

■ 

324 

-49.4 

1.65 

1.17 

1.0 

mmm 

NON STA. 

48 



730906 

211947 

212000 

IjD 

100 

-4.71 

1,45 

.318 


.211 

16.19 





212003 

212134 

H 

650 

-7.58 

.931 

8.57* 

.808 

0-.060 

NON STA. 



5 OAS-3 


212142 

212236 

■ 

352 

-10.8 

1.24 

3.21 

1.0 

HBfBB 

NON. STA. 









































































































































TABLE 3.8 <CONT. } 


WAP 

PASS 

MODE 

DATE 

START 

STOP 

PbiNT 

<OcGS> 

SAW*. 

SIZE 

»EAN 

(MTRS) 

SIGKA 

<MTRS) 

F 

STATISTIC 

NORM. 

TEST 

SIS. 

FREQ. 

BA.\OS 

(KZ> 

SIGNIFICANT 

AOTOCOPREUTION 

COEFFICIENTS 

SO 

28 

S OAS-2 

730907 

203582 

203712 

<.6* 

720 

-2.61 

.828 

.306 

.731 

0-.065 

NON STA. 



5 DAS-3 


203715 

203805 


326 

-3.99 

.876 

.161 

.535 

0-.121 

NON STA. 



5 OAS-1 


203516 

203529 


100 

-.756 

1.22 

.165 

.669 

.567 

13.16,16 

92 

57 

5(nOAS-2 

731202 

182632 

182606 

<.5* 

650 

38.0 

3.67 

..79* 

1.0 

0-.600 

NON STA. 

93 

57 



182705 

162715 

<.S* 

76 

19. 1 

1.78 

.186 

.693 

NONE 

NONE 


71 

5<*>0AS-1 

7^f^:oi 

133237 

133251 

.5* 

106 

17.56 

1.72 

.720 

.5X6 

.072.. 176 

1.2.5.11.13 

1 


5(6)OAS-2 


133256 

133633 


706 

20.6 

3.B9 

6,65* 

1.0 

1 

0-.067 

i;ON STA. 

132 

1 

5<2)OAS-*l 

760109 

156969 

156959 

■ 

76 

57.7 

1.36 

.175 

.950 

- ^ 

.622 

3.6.9.10 

13.16.20 


1 

5<2)0AS-2 


155002 

155161 

■ 

706 

67.1 

3.52 

5.86* 

■ 

0-.067 

NON STA. 

S33 

86 

S OAS-I 

760120 




m 

66.3 

1.63 

1.26 

.560 

.659 

189-12.20 



5 OAS-*2 


191830 

192010 

1 

720 

82.8 

1.56 

1.06 

■ 

0-.217 
.306-. 369 
.656-. 521 
.760-. 825 

1-6. 6-9 










































































TABLE 3.9 


MODE 5 MODERATE POINTING ANGLE 



r 


1 

PASS 

MODE 

DATE 

START 

STOP 

POINT 

COEGS) 

SAMP. 

SIZE 

MEAN 

<KTRS) 

SIGVA 

(MTRS) 

F 

STATISTIC 

■Si 

SI6. 

FREQ. 

BANDS 

<HZ) 

SIGNIFICANT 

AUT0C0FRELATI0T4 

COEFFICIENTS 


B 

5 OAS-1 

730609 

151401 

151546.5 

mi 

114 

91.3 

1.50 

.70 

.656 

NONE 

NONE . 

18 

11 

S(2)OAS*l 

730614 


145415 

.7* 

112 

60.8 

2,25 

2.54 

.975 

.1268.153 

.767 

182.3.15.16 

19 



5(2)0AS*2 


145417 

145516 


352 

62.8 

1.72 

■ 

1 

0-.178 

NON STA. I 

37 

21 

S(3)0AS-1 

730901 

153022 

153036 

.75* 

108 

-38.7 

1.47 

1.24 


NONE 

1.4. 6. 7 

73 

39 

5( 1 )0AS-1 


194549 

194602 


100 

-54.6 

2.31 

.760 

.762 

.234 

i-4.9.13 



SC 1 lOAS-3 

730913 

194747 

1.94836 

. 88 * 

347 

-84,1 

■ 

13.9* 

1.0 

0-.114 

NON STA. 

74 

39 

5(2)0AS-1 

730913 

194927 

194940 


100 

-93.2 

1.85 

.648 

.606 

.148 

8.20 



5C2)0AS-2 


194946 

195122 

. 88 * 

686 

-89.8 

2.47 

3.11 

1 

0-.057 

NON STA. 


1 

5<2)0AS*3 


195125 

195216 


342 

-100.3 

1.83 

24.5* 

.999 

0-.114 

NON STA. 

105 

65 

5 OAS-1 

731215 

000100 

000113 

.75* 

100 

-66.7 

2.10 

1.12 

.576 

NONE 

1.4 

106 

67 

5 OAS-1 

731218 

020019 

020030 

.75* 

84 

36.1 

n.94 

.414 

.733 

1.02 

18384.12 

































































































































TABLE 3. 9 (CONT. ) 


mp 


MODE 

DATE 

STpm 

STOP 

POINT 

(DECS) 

■ 

MEAN 

(MTRS) 

SIGMA 

(MTRS) 

F 

STATISTIC 

■n 

SIG. 

FREQ. 

BANOS 

(HZ) 

SIGNIFICANT 

AinrOCORRELATION 

COEFFICIENTS 

loe 

1 

S OAS-2 


020033 

020193 


966 

42. S 

2.60 

2.10 

1.0 

0-.113 
.207-. 309 
.579-. 627 

NCN OTA. 

109 

1 

S<U0AS-2 

731210 

119011 

119690 

.TO* 

1 

22.0 

3.90 

1.01 

.978 

0-.057 

NON STA. 

110 

60 

9<2)0AS-1 

731210 

119013 

119627 

.70* 

109 

-.199 

2.09 

1.17 

.099 

1.20 

2 

S(2)0AS-Z 


114«30 

110010 


600 

7.21 

2.17 

9.02* 

• 996 

0-.060 
.216-.399 
.933-. 060 
.699-. 761 
.660-. 925 

NON STA. 

113 

71 

Sf2)DAS-l 

790101 

132017 

132032 

.70* 

119 

30.91 

1.71 

1.01 

.782 

.911 

1.16.20 

S(2)0AS>2 


132039 

132719 


720 

4S.4 

1 

1.63 

9.03* . 

1.0 

0-.217 
.360-. 990 
.606-. 790 
.636-. 901 

i-4.e-i2, 

16-20 

119 

Tl 

S(3)0*S>1 

790101 

132009 

132913 

.70* 

109 

33.0 

1.93 

.731 

.966 

.901 

1.0.16 

9yt3)OA$-2 


132916 

133009 


702 

. . .. 

43.2 

1.00 

2.03 

• 091 

0-.067 

NON STA. 

120 

70 

9<t>tMS-l 

790100 

16340S 

163916 


69 

06.2 

1.91 

2.39 

.776 

NONE 

1 

SOIOAS-2 

163920 

1 

163000 

.60* 

709 

09.7 

3.30 

2.06 

1.0 

0-.199 

NON STA. 





























table 3.9 (CONT.) 


u» 

I 


vO 


*AP 

PASS 

KCCE 

DATE 

START 

STOP 

POINT 

(OEGS) 

SAWP, 

SIZE 

tCAN 

(t^RS) 

siGm 

(MTRSt 

F 

STATISTIC 

NORM. 

TEST 

SIG. 

FREQ. 

BANOS 

(HZl 

SIGNIFICANT 

AUTOCORRELATION 

COEFFICIENTS 

129 

78 

S(2)OAS-2 

740108 

164444 

164623 

.75* 

704 

-33.3 

6 



2.01 

.998 

0-.0655 

NON STA. 

I^S 

63 

5(UDAS-1 

740114 

153433 

153452 

.75* 

1C4 

2.03 

2.21 

.748 

.557 

.826; 1.20 

1,«.19.20 

146 

83 

S(2)DAS-2 

740114 

153841 

154021 

.75* 

704 

-20.7 

2.72 

9. 78* 

.965 

0-.067 
.388-. 610 

NON STA. 

D 

8S 

5 OAS-l 

740118 

204835 

204848 

.75* 

100 

78.4 

2.02 

1.45 

.734 

.391 

1.15 

216 

1 

97 

5(1)045-1 

740131 1 

155628 

155641 

.75* 

92 

-5.14 

2.40 

.986 

.944 

.313 

1.9.10.16 



SCl)DAS-2 


155643 

155823 


650 

8.37 

2.55 

3.19 

■ 

0-.421 

NON STA. 

217 

97 

S(2>0AS~1 


155993 

155956 

.75* 

too 

12.6 

1.69 

.387 

.620 

1.56 

ltlS.18 



S(2>0AS-2 

■ 

155959 

160137 


650 

33.3 

3.75 

. t.l2 

.989 

0*-.132 

NON STA. 
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TABLE 3.10 


MODE 5 HICH POINTING ANGLE 


1 

PASS 

MODE 

DATE 

START 

STOP 

POINT 

(DcGS) 

SAMP. 

SIZE 

KEAN 

tMTRS) 

SIGMA 

{MTRS) 

F 

STATISTIC 

NORM. 

TEST 

SIC. 

FREQ. 

BANOS 

(H2) 

SIGNIFICANT 

autocorrelation 

COEFFICIENTS 

34 

15 

5 OAS-2 

730812 

024657 

024829 

l.l* 

600 

-23.5 

3.56 

9.31* 

.530 

0-.534 

1-6.8-16 

19.20 

54 

32 

5n)DAS-l 

730901 

130616 

13C628 


Bi 

23.9 

8.16 

14.3* 

.544 

.170-. 425 

1-5.8-17 



5( t )0A3-2 

. . 


130639 

13C810 

1.6* 

768 

19.3 


20.9* 

1. 

0-.345 

1-12.15-20 



5{ 1 )0AS-3 


130815 

1 30902 


360 

12.8 

9.33 

27.3* 

.997 

0-.174 

1-14.18-20 

55 


5(2)OAS-l 

7309U 

130946 

130957 

D 

84 

4.35 

12.65 

42.9* 

wm 

.279-. 372 

1-5.9-19 



5(2)DAS-2 


131001 

131114 

■ 

756 

-7.05 

6.64 

49.9* 

m 

0-.279 

1-12 

16-20 

■ 


S(2)0A5-3 


131143 

131232 

■ 

360 

-6.13 

5.25 

4.43* 

.993 

0-.297 

1-1 1 
16-20 

Q 

35 

5( 1 )OAS-2 

7309X2 

125544 

122727 

1.4* 

759 

-5.04 

6.42 

22. !• 

B9 

1 

0-.335 

i-12, 16-20 

■ 


5U)dAS-3 


122730 

122620 


384 

rai 

4.87 

3.14 

BI 

0-.346 

1-10.16-20 

■ 

i 

5(1)0AS-1 

730912 

1 

165215 

165228 


92 

50.0 

7.99 

3.19 

.999 

.255-. 340 
.510 

1-4,7-XT 



5(nOAS-2 


165231 

165411 

1.25* 

650 

52.3 

3.40 

4.94* 

.981 

0-.176 
.361-. 565 

all 

1 


■ 


165413 

165503 


200 

i 

58.5 

3.90 

1.53 

.623 

0-.0B37 
.167-. 251 
.335-. 530 
.614-. 642 

1-4.8-10. 

15-18 

ll 


S(2)0A$-I 

730912 

165553 

165559 

m 

46 

62.1 

6.99 

3.50 

.996 

.340-. 679 

■SI 
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TABLE 3,10(CCNT.) 


! 



DATE 

START 

STOP 

POINT 

(OCGS) 

SAM*. 

SIZE 

MEAN 

(MTRS) 

SIGMA 

(MTRS) 

F 

STATISTIC 

1 

SIG. 

FREQ. 

BANOS 

(HZ) 

SIGNIFICANT 

AUTOCORRELATION 

COEFFICIENTS 

62 

36 

S(«)CM»>2 

730912 

t6S«09 

165737 

1.3* 

578 

63.2 

3.36 

2.95 

.997 

0-.392 
.541*-. 595 

1-7.9-16 



5(2lOAS-3 


I9S7S1 

165835 


240 

65.5 

3.0* 

2.57 

.996 

.0651-. 293 
.391-. 423 

1-7.12-20 

63 


S(3>CMS-1 

730912 

>65919 

165930 

H 

84 

62.4 

6.69 

5.79* 

• 939 

.093-. 186 
.372-. 465 

1-4.9-15 


■ 

5(3)0AS-2 


165932 

170112 

■ 

650 

53.4 

5.60 

6.20* 

D 

0-.276 

1-12.15-20 

B 

25/36 

S(«)OAS-2 

730912 

171244 

171423 

■B 

650 

46.6 

8.82 

16.4* 

m 

0-.348 

NON STA. 

66 

37 

5 DAS-1 

730912 


201049 


110 

70.4 

11,94 

2.54 

.969 

0-.213 

.855 

MJL 



5 DA5*2 


201056 

201232 

1.4* 

640 

76.8 

5.81 

23.2* 

D 

0-.354 

1-14,18-20 



5 OAS-3 


201234 

201325 


294 

84.6 

4.46 

10.6^ 

DSi 

0-.239 

AU. 

B 

» 

5(2)065-1 

730913 

160952 

181003 

1,55* 

84 

41.6 

9,27 

10.53* 

.995 

.093,. 279 


99 

61 

5 OAS-1 

731205 

162108 

162120 

l.l* 

92 

-683 

6.17 

7.86* 

.995 

.340,. 510 

1-3.6-12 



5 D6S-2 

i 

162120 

162222 


468 ' 

2.71 

6.91 

13.25* 

■ 

0-.05 
.133-. 384 
.534^.564 


136 

79 

S(MOMS-l 

7M>109 

160044 



too 

-68.1 

1.93 

1.86 

.891 

NONE 

1,16 



5(4)065-2 


160100 


1 

650 

-54.3 

2.34 

2.89 

.916 

0-.060 
.144-. 204 
.288-. 349 
.505-. 565 
.793-. 925 

NON STA. 













































































































































within thn progr 
arc in eoi^tttnd ns 


where {X^} ere the altioeter range residuals, X is the residual s«an and H is 
the saaiple slse. Thus o is a neasure of the dispersion of the data about the 
■can. Table 3.11 displays the overall residual slgnas for each ■ode/polnting 
«n gi<! range. Because of the snail sasqile slses data for Mode-3 low and aoderate 
pointing angles has been grouped together. Figures 3.1, 3.2 and 3.3 show the 
distribution of residual slgnas in nodes 1, 3 and 5 for each pointing angle 
range. In each case the overall signs for each siode/pointing angle is defined 
as 


SAP tha residual standard deviation, for a given data 


r 2 „=2 

I X. - NX 
fl ^ 


iI-1 


1/2 


(3.1) 


6 




(3.2) 


where {o^} are the arc slgnas and {M^} are the arc sanple sizes. In coi^ting 
6 for each dau set obviously outlying were suppressed. 


Figure 3.4 is a scatter diagram of residual signs as defined in equation 
(3.1) versus tracker determined pointing angle. For pointing angles above .9 
degrees each individual data point is shown. Due to tha large number of data 
points below .9 degrees only tha range of residual slgnas, with obvious out- 
lyars deleted, is shown. A least squares regression line was fitted to tha 
52 data points with pointing angles greater than .9 degrees yielding tha 
statistical relationship t 

Residual Sigma - 3.38 meters 7. 34* (Pointing angle in degress) 

The corresponding correlation coefficient is .5306. 
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TABLE 3.11 

OVERALL RESIDUAL oioMAS 
(METERS) 


POINTING 

ANGLE 

RANGE 


MODE- 3 

mode-5 

LCM 

1.S6 

1.64* 

1.27 

MODERATE 

2.43 

2.18 

HIGH 

5.51 

5.28 

5.68 


* INCLUDES BOTH LOW AND MODERATE 
POINTING ANGLE DATA. 
















FIGURE 3.1 


OISmiBUTION OF MODE 1 
RESIOUM. STANDARD DEVIATION 
IN METERS 
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012 34 S6789 

LOW AND MODERATE POINTING ANGLES 


jO * 5.28 



HIGH POINTING ANGLES 






RESIDUAL SIGMA METERS 


FIGURE 3.4 
RESIDUAL SIGMA 
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jOC 


Figure 3.5 Is a typical plot of the grouped power spectral coefficients as 
described in Section 2.2.1. The frequency tabulated is the group center fre- 
quency. For example the frequency band associated with the first grouped co- 
efficient (.0333 Hz center frequency) is zero Hz to .666 Hz. In this example 
the grouped coefficients had 8 degrees of freedom. The 99th percentile of the 
chi square distribution with 8 degrees of freedom is 20.1. The first two grouped 
coefficients exceed this value, therefore the significant frequency band is 0 Hz 
to .133 Hz. In this example there is a well defined low frequency signal and 
the measurement noise power is close to being xinlformly distributed over the 
frequencies above .300 Hz. 

Figures 3.6 through 3.13 are bar charts showing the location of significant 
frequency bands for each mode/polntlng angle range. As discussed in reference 
[1], recoverable geoid features are contained in that part of the frequency 
spectrum below .164 Hz. In spite of having selected data arcs without prominent 
geoid features displayed on the residual time series plots, figures 3.6 through 
3.13 indicate that most of the data arcs processed did contain significant power 
spectral components below .164 Hz. For the majority of arcs having low or 
moderate pointing angles these bands are well separated from significant higher 
frequency bands. This indicates that signal (i.e. geoid) and noise ccHoponents 
of the range residual data can be separated with appropriate filtering techniques. 

For high pointing angle arcs (figures 3.8, 3.10 and 3.13) this separability 
is not present. The measurement noise power is concentrated at lower frequencies 
which tend to overlap the 0-.164 Hz band containing recoverable geoid features. 
From reference [1] a reasonable explanation for this behavior is that at higher 
pointing angles the received waveform changes shape causing the altimeter 
tracking loop which acts as a low pass filter to become more "slugggish". 
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FIGURE 3.6 
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FIGURE 3.8 

DISTRIBUTION OF SIGNIFICANT 
FREQUENCY BANDS 
NODE 1 HIGH POINTING ANGLES 


• MAP 


MAP 

MAP 

MAP 


6 

DAS 


DAS 


DAS 

28 

DAS 


DAS 

33 

DAS 


DAS 

108 DAS 


DAS 



l.QHZ 


1.5HZ 


2.0HZ 


MAP too DAS 1 
MAP U8 DAS I 



1 1 1 


FIGURE 3.9 
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FIGURE 3.12 
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3.3 R»f«r»ac>« 

1. McGoogan, Laltao and WallSt Suaaary of SKYLAB 8»193 Altlaatar Altltnda 
Results. NASA TM X-69355, NASA Wallops Flight Center, February 1975. 
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In th« previous section (3.0) » sunaaery statistics have been presented for 
approxinately 184 arcs of S-193 altiaeter residuals. For a typical data arc 
these residuals can be modeled as 


r(t) - s(t) + n(t) 

Where s(t) is a low frequency, non-statlonary signal due primarily to geoid 
variations and n(t) is stationary, broadband aMasuremsnt noise. This section 
describes detailed evaluations of the statistical properties of the noise pro- 
cess, n(t), for twenty-seven data arcs (Table 4.1) which were selected to span 
all data acquisition altimeter modes/submodes and the range of pointing angles 
(low, moderate, high). This evaluation was done by 1) selecting data arcs where 
s(t) was (Avlously small and/or at low frequencies, 2) in some cases high pass 
filtering the data to further attenuate s(t) and 3) coq>uting deacrlptive ata- 
tlstlcs within the RAP program. The high pass filter, used to remove the signal 
s(t) , was obtained by low pass filtering the data with a 99 point quadratic 
midpoint filter and taking the high pass output as the difference between the 
filter midpoint value and the corresponding input data point. CALOCMP plots of 
the resultant time series, power spectrum end autocorrelation function were 
obtained using the optional RAP plot tape output. Also a 20th order autore- 
gressive model of the noise process was formulated and the corresponding analytical 
power transfar function plotted. 


In Section 3.2 it was shown that, for high minting angles off nadir (i.e. 
greater than .9*) the frequency bands occupied by tha signal s(t) and the noisa 
n(t) show considerable overlap. Examination of the unfiltered data time series 
plots for all of the high pointing angle data arcs in Table 4.1 (Maps 6, 12 
and 66) showed that s(t) for Map 6 DAS-2 Includes the largest trend which is 
approximately .22 meters per second (Figure 4.1.29A). With the trend removed 
by a 99 point high pass filter the data standard deviation is decreased from 
6.23 meters (unfiltered) to 5.30 meters (filtered) (Figure 4.1.29B). The power 
spectrum of the unflltered data (Figure 4.1.30A) shows that the spectrel com- 
ponents due to this trend (0 to .10 Hs) are smaller by approximately 4db than 
the dominating components due to the noise oscillations induced within the 
tracking loop (.20 to .26 Hs). Comparison with the filtered date power spectrum 
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TABLE A.l 


RESIDUAL DATA ARCS SELECTED 
PGR DETAIL ANALYSIS 


mp 

PASS 

MODE 

SUBMOOE 

TRACICR KTERMltCO 
POINTING 

39 

1 

1 

OAS-1 

DAS-2 

DAS-3 

<.6* (LOW) 

205 

97 

1 

DAS-1 

DAS-2 

DAS-3 

.7S* (MODERATE) 

6 

7 

1 

■33m 

1.2* (HIGH) 

1 

24 

3 

■Sm 

<.6* (LOW) 

10 

8 

3 

DAS-l 

DAS-2 

DAS-3 

.75* (MODERATE) 

12 

9 

3 

DAS-1 

DAS-2 

DAS-3 

1.1* (HIGH) 

47 

25 

5 

DAS-1 

DAS-2 

DAS-3 

<.S* (LOW) 

7A 

39 

5 

mmm 

B3m 

.85* (MODERATE) 

M 

37 

5 

DAS-1 

DAS-2 

DAS-3 

1.4* (HIGH) 

1 
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(Figure 4.1.30B) showed that at fraquancles above .25 Hs the ahape of the 
filtered and unflltered power apectra are Identical. Thus the presence of this 
trend had only a 15X effect on the data standard deviation and did not doadnate 
the power spectral plot. For these reasons It was decided to present tine series 
plots, power spectral, autocorrelation functions and autoregressive transfer 
functions derived from unflltered data for the nine arcs with large pointing 
angles off nadir. (Haps 6, 12 and 66). In addition, the 99 point filter wes 
not applied for data arcs shorter than 30 seconds. This Included Mode 3 DAS-1, 
Mode 3 DAS-2 and Mode 5 DAS-1. 



4.1 Mode 1 Results 


Ftob Table 2.1 the three mode 1 data acqulaltlon aubmodea (D4S-1, DAS-2, 
and DAS-3) are 52 aeconda* 63.44 aeconda, and 63.44 aeconda long reapectlvely. 
During DAS-3 the antenna pointing la offaet .431* in pitch to evaluate the 
effect of a known change in pointing angle. Table 4.2 diaplaya aumaary ata- 
tiatica for the 9 aelected mode 1 data area. The correaponding time aerlea 
plota» power apectra» autocorrelation functlona and autoregreaalve model power 
tranafer functlona are diaplayed in Figurea 4.1.1 through 4.1.36. 

Aa would be expected with increaalng pointing angle the nolae atandard 
devl<itlona (Table 4.1 - Filtered Data) Increaae and the power apectra ahow a 
well defined attenuation at the higher frequenclea. 

At low pointing anglea (Map 39) the meaaurement nolae » n(t), reaend>lea a 
white noiae proceaa for all 3 aubmodea (DAS-1, DAS- 2 and DAS-3). The power 
ap.;ctra (Figurea 4.1.2» 4.1.6 and 4.1.12) are relatively flat and the correa- 
pondlng autocorrelation functlona (Figurea 4.1.3, 4.1.7 and 4.1.11) have few 
coefflclenta which are algnlficantly different from zero. The atatiatical 
characterlatlca of the data ahow no apparent change due to the .431* pitch off- 
aet during DAS-3. 

At moderate pointing anglea (Map 205) during aubmode DAS-1 the power 
apectrum (Figure 4.1.14) and the autocorrelation function (Figure 4.1.15) 
indicate little or no departure from a white noise proceaa. Submodes DAS-2 
and DAS-3 show significant departures from a white nolae process. The power 
apectra (F:lgurea 4.1.18 and 4.1.22) show considerable attenuation at the higher 
frequencies with the DAS-2 spectrum starting to taper off between 1.5 and 1.6 Hz 
and the DAS-3 spectrum starting to taper off between 1.25 and 1.35 Hz. This 
difference in break frequencies could possibly be attributed to the .431* pitch 
of .set of the altimeter antenna during DAS-3. Comparison of the autocorrelation 
functions for DAS-2 and DAS-3 (Figures 4.1.19 and 4.1.23) shows that the DAS-3 
data has more autocorrelation coefficients significantly greater than zero. 

At high pointing angles (Map 6) the time series plots for all 3 submodes 
(Figures 4.1.25, 4.1.29 and 4.1.33) indicate that the noise process is dominated 
by low frequency oscillations which are approximately 23 meters peak to peak 
during DAS-1 and DAS-2 and as high as 35 meters peak to peak during DAS-3. The 
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power spectral plots (Figures 4.1.26, 4.1. 30A end 4.1.34) show that the prin- 
cipal oscillation frequency for DAS-1 is at .4 Hs, for DAS-2 between .19 Hi and 
.25 Hi and for DAS-3 between .22 and .35 Hi. The .431* pitch offset during 
OAS-3 causes the noise sigma to increase from the DAS-2 value of 5.3 meters 
(Table 4.1) to 8.04 meters. Cooq;>arison of the DAS-2 and DAS-3 power spectral 
plots (Figures 4.1.30A and 4.1.34) show a much higher concentration of energy 
at the lower frequencies for DAS-3. At 1.0 Hs the DAS-2 power spectrum (Figure 
4.1.30A) has decreased about 5 db from its peak value while the DAS-3 power 
spectrum (Figure 4.1.34) has decreased about 18 db. 


TABLE 4.2 


MODE 1 DETAIL ANALYSIS 


mp 

PASS 

MODE 

POINT 

(OEGS) 

SAMP. 

SIZE 

MEAN 

(MTRS) 

SIGMA 

(MTRS) 

F 

STATISTIC 

39 

22 

*1 DAS-l 


208 

17.9 

1.16 

1.28 



*1 OAS-2 

<.6* 

408 

22.5 

1.03 

1.36 



•l DAS-3 


408 

23.3 

1.12 

.362 

205 

97 

*1(6)DAS-1 


280 

0 
. 

CD 

1 

1.89 

1.35 



1(6)DAS-2 

.75* 

308 

-46.1 

1.69 

.583 



1(6)0AS-3 


368 

-45.4 

1.66 

4.78 

6 

7 

1 DAS-1 


192 

89.0 

5.30 

1.16 

h ■ ■ ■ ■ 



1 OAS-2 

1.2* 

368 

94.5 

6.23 

1.52 



1 DAS-3 


352 

146.3 

8.55 

10.13 


♦indicates data filtered by 99 POINT QUADRATIC MIDPOINT HIGH 


NORM 

TEST 


SIGNIFICANT 

AUTOCORRELATION 

COEFFICIENTS 


.446 


.929 


.995 


.960 


.846 


.082 


.796 


.101 


.995 


SIG. 
FREQ. 
. BANOS 
(HZ) 


.676--. 714 


2.37-2.43 


1.53-1.59 


NONE 


.228-. 355 
.457-. 583 
.913-. 964 
1.37-1.42 


.191-. 223 
.382-. 488 
.764-. 807 
.891-. 934 


0-.081 
.407,. 529 
.651-. 692 
.895 


0-.488 


0-.510 


3,5 


8,12,16 


8 


NONE 


1,9,10,14 


1,2,11-14,19 


1,2,15,16 


ALL 


ALL 


ss filter 


I 

i 
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FIGURE 4.X.1 
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FIGURE 4.1.3 


MAP- 39 PASS 22 MODE 1 OAS-1 
AUTOCORRELATION FUNCTION 
POINTING a .6* 
(FILT ER ED) 








FIGURE 4.1.4 

AUTOREGRESSIVE MODEL 

MAP- 39 PASS-22 MODE- I DAS- I 
(FILTERED) 


GAIN 

FREQUENCY (FRACTIONAL) 


PLOT OF 

POVfER TRANSFER FUNCTION 08 


U, C.7C?‘0705('t:-01 

0,766E-01 C.(‘j197i04E-01 

0.153E 00 0.li>‘>57567 

O.rSOF 00 0.t>^7l0391 

0.306E 00 " 2,9166044 “ 

0.38JE 00 0.046957^4 

0.460E 00 Cl. 47879834 

0.536F 00 ' 0.57690419 

0.6136 00 1,4??4263 

0,6096 00 3,0171069 

0.7666 00 ■ 1,400778b 

0,0436 00 0.79065071 

0.9196 00 0.(7077034 

0,9966 00 1,6055355 

0.107E 01 2,47953.'>6 

0,1156 01 1,3793145 

0.1236 01 0.91731311 

0.130E 01 1,0036881 

0.1386 01 1.6607148 

0,1466 01 2,12(6672 

0.1536 01 1,203431? 

O.lMf 01 0,77794597 

0,1696 01 0.74672476 

0.1766 01 0.99710283 

0,1846 01 1.3,'94V96 

0,l9lE 01 1 ,?c39295 

0.1996 01 1,0069308 

0.2C76 01 0.99353943 

0,2146 01 1,2974(74 

0,2226 01 1,7425622 

0,2306 01 1,5564592 

0.2376 01 1.1457564 

0,2456 01 0.96976595 

0,2536 01 1,0164174 _ 

0,2606 01 1.C595040 

0,2M6 01 (.99(54124 

0.2766 01 0.(951535? 

0,2836 01 0.i'254?043 

0.2916 01 0.(1519644 

0.29:6 01 0. <9001619 

O.SC’.C 01 0.57332317 

0,314? 0.’. 0.52500668 

0.322F 01 C.5BH3?O02 ' 

C.329F 01 0.(1)02658 

0.337F 01 1, 2057633 

0,3<:E 01 1,4648067 

0,3526 01 1.4530082 

0.3M.E 01 1,5501126 

0.36(6 01 1,6645970 ' 

0.375E 01 1,9365965 

C,303E 01 1,5389435 


0 



20 
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FIGURE 4.1.5 

M4P-39 PASS 22 MODE 1 DAS-2 
ALTIMETER RANGE RESIDUALS 
POINTING s .6* 
(FILTERED) 
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FIGURE 4.1.6 

MAP-39 PASS 22 MODE 1 DAS-2 
POHER SPECTRAL DENSITY 
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FIGURE 4.1.7 

MAP-39 PASS 22 MCX3E 1 DAS-2 
AUTOCGRRBJkTION FUNCTION 
POINTING s .6* 
(FILTERED) 




FIGURE 4.1.8 

AimiREGRESSIVE MODEL 


MAP-39 PASS-22 MODEL-1 OAS-2 
(FILTERED) 


GAIN 

FRBOUENCY (FRACTIONAL) 


PLOT OF 

POMER TRANSFER FUNCTKM OB 


0. 


0.829923766-01 

0.764E- 

-9i 

1^.962030386-01 

0.153E 

9u 

0.13359478 

0.23JE 

SU 

0.36**12264 

0.394c 


1,5404295 

0.383E 

90 

3,9245282 

0.44.SS 

9U 

1,424»0.>9 

0 .S36c 

9*1 

1.0702/51 

0.613E 

9u 

1,0.«92504 

0.639E 

90 

0.94/50017 

0.764E 

00 

C.87537881 

0.6<3E 

90 

1,0270576 

o.v:9t 

90 

1, 7154722 

0.9956 

?i) 

2,4453704 

0.107K 

9i 

1.2707627 

O.llSfc 

9X 

C.e»092/fi7 

0.123I; 

91 

0.54844112 

0.l3vE 

n 

C.65923M6 

0.1336 

91 

1,1005585 

0.l<3e 

91 

2,:*5379.)2 

0.tS3E 

91 

1,9855592 

O.lOlE 

91 

l,2779<?7 

0.1696 

)1 

1.06225V6 

0.17t.F 

01 

1.0559799 

0.1846 

n 

1.0RC1566 

0.1916 

31 

C.K4666204 

0.l99h 

31 

(,./6'»34/86 

0,2076 

91 

0.oC775i9.I 

0.2146 

91 

0.59549073 

0.2226 

91 

1.2341CH4 

0 .?3i)E 

91 

1,3781604 

0.2376 

01 

1,5286223 

0.2456 

91 

1,0255263 

0.2536 

91 

1.79R50O1 

0.2606 

01 

1.16995:34 

0.2606 

n 

C,.7«289762 

•).2745 

n 

0.«. 724056 7 

0.2836 

01 

1). 11037423 

0.2916 

Cl 

1,3754431 

0.2996 

91 

2.2954448 

0.3946 

01 

l.f 186290 

0.314f 

01 

t..», <748435 

0.322? 

01 

p. 75804694 

0.3296 


i..» 2356882 

0.3376 

?1 

I.P0157.36 

0.3456 

r; 

1,3382158 

0.3526 

Cl 

1.2715444 

0.3A.6 

03 

l.(8098O« 

0.34*.-6 

91 

i.V7S07546 

0.3756 

91 

r.V.»721264 

0.303% 

Cl 

O.t 1463404 
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FlUURE 4.1.9 

MAP-39 PASS 22 MODE 1 DAS-3 
ALTIMETER RANGE RESIDUALS 
POINTING = .6* 
(FILTERED) 



20.0 
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FIGURE 4.1.11 




FIGURE 4.1.12 

AUTOREGRESSIVE MODEL 

MAP- 39 PASS-22 MODE-1 DAS- 3 
(FILTERED) 


GAIN 

FREQUENCY (FRACTIONAL) 


0 , 

O,766E«»01- 
0.153E 00 
0i230E 00 
0.306E 00 
0.383F. 00 
0,460F 00 
0,536L 00 
0,613r- 00 
0,689E 00 
0.766E 00 
0.R43E 00 
0,919E 00 
0,996E 00 
0,107E 01 
0.U5E 01 
0,123E- 01 
0,130E 01 
Q , 1 3 8 E 01 
0,1-^6E 01 
0,1?>3E 01 
0il6lE 01 
0,1£>9E 01 
0.176F- 01 
0,1 84E 01 
0.191E 01 
0.199E 01 
0,207E 01 
0,2140 01 
0,22?E 01 
0.230E 01 
0.237E 01 
0,245t: 01 
0,2^3E 01 
0,?£>0E 01 
0,268E 01 
0,276E 01 
0,2a3fe 01 
0.291E 01 
0,299E 01 
0.306E 01 
0.314E 01 
0.322E 01 
0,329E 01 
0,337R 01 
0.345E 01 
0,352E 01 
0,360E 01 
0,368E 01 
0*37!1F 01 
0,3 o3l 01 


0.64973102E-01 

O,7039642eE-Oi 

0.1299761B 

0.34717964 

1,7776574 

1,564R26& 

0.69^75690 

0.63600672 

0,95352429 

1,5589097 

l,2V?5901 

0.6981 9*4 6 6 
0,66025806 
1,1022711 
1,3455264 
1 ,11..0246 
0.67667951 ' ' 

0,90i;53ao7 
1,3102206 
2,3203524 
2,6459874 
1,8210084 
1 ,4060159 
1,2991104 
1,2028545 
1,0236745 
0.90946935 
0.95997937 
1,1979052 
1,4284756 
1,2692974 
1,0126642 
0.94552538 
1,0861984 
1,3195536 • 

1,2659644 
0,98709608 
0,63075614 
0,65215970 
1,0337843 
1,2683063 
1,309094V 
1,2075378 
1,3.980255 
1,3540298 
1,5634918 
1,5101753 
1 ,2151519 
0.97596747 
0.06296541 
(1,83483116 
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FIGURE 4.1.13 


MAP-205 PASS 97 MODE 1 OAS-1 
ALTIMETER RANGE RESIDUALS 
POINTING = .75® 
(FILTERED) 












FIGURE 4.1.15 

MAP-205 PASS 97 MODE 1 OAS-1 
AUTOCORRELATION FUNCTION 
POINTING = .75* 


(FILTERED) 





FIGURE 4.1.16 

AUTOREGRESSIVE MODEL 


FREQUENCY 


• GAIN 
IFRACTIONAL) 


MAP-205 PASS-97 MOOE-1 OAS-1 
(FILTERS)) 

PLOT OF 

POWER TRANSFER FUNCTION OB 


i 


0 . 


0.766E- 

■01 

0.1 53E 

00 

0.?3oC; 

<)0 

O.3O0E 

00 

0.303K 

•)0 

0»46iiE 

00 

0.53/.F. 

00 

0.613E 

00 

0.6H9E 

00 

0.766F 

00 

0.e43E 

00 

0.919E 

00 

0.996E 

00 

0.li)7r: 

0 ; 

O.UfE 

ftj 

0.123E 

•01 

0,13‘j£ 

01 

C.136F 

9* 

0.1 46H 

01 

0.153E 

01 

O.lcii 

0* 

C.169E 

01 

0.176E 

0 ; 

0,1«4E 

01 

O.lSlS 

r * 

0.179c 

01 

0.207E 

01 

0.214c 

0 : 

0.22PE 

ni 

0.230E 

31 

0.237? 

01 

0.?46F 

01 

0.253F; 

cl 

0.26r.f: 

31 

0.?6P(: 

Oi 

0.27*-t 

n 

0,2r.3E 

01 

0.?9lE 

Cl 

0.?59r; 

01 

0.3f6h 

01 

0.3J4E 

31 

C.322r 

31 

0.329E 

31 

C.337F 

01 

0.346? 

ri 

C.35LE 

Ol 

0.36(‘fc 

01 

0.3(6? 

01 

0.3750 

f J 

0.3».3t 

0 : 


0.*73<1<2& 

0.i'fU7i79 

t.7J3f3897 

i,6D90^74 

1 ,\»<7c36 
0.^^437420 
l.l«?1047 

1. V(i9«p{»l 

2, «‘*o2577 
1.3<>A38ti7 

1, K(67?a 
l,3?C-2949 
1,5349575 
1,2319)63 

t.V2646;j53 

Ii.V3962^»*<5 

1.4540350 

2, J>4?f,7<.C 
2.0fe'G66.i7 
1,0034790 

( .72f 34748 
C.7/232-.54 
1,1086376 
1,4791297 
1, '<769/37 
()./3{>74 686 
0 ,o2070o35 
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C. >9907255 
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) ,1.327539 
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FIGURE 4.1.17 

MAP-aOS PASS 97 MODE 1 OAS-2 
ALTIMETER RANGE RESIDUALS 
POINTING s .75* 
(FILTBED) 





20.0 











FIGURE 4.1.19 

MAP-205 PASS 97 MODE 1 DAS^2 
AUTOCORRELATION FUNCTION 
POINTING = .75® 


(FILTERED) 





FIGURE 4.1.20A 


AUTOREGRESSIVE HQOEl. 

MAP-205 PASS-97 MOOE-1 OAS-2 
(lAFILTERED) 


GAIN 


PLOT OF 


FREQUENCY 

(^ACTIONAL) 

0 • 


4'.0444299 

0t766E< 

-01 

3.4119005 

P.153E 

00 

2'.7346469 

0.230E 

00 

2.6085774 

0.306E 

00 

3.8937633 

0.383E 

00 

5.5856865 

0.460E 

00 

4.5402775 

0.536F 

00 

2.8416538 

0.6t3E 

00 

2:i79c8<8 

D,«a9E 

00 

2.136H77J " 

0,7<)6E 

ou 

2.4C8V710 

P.e43E 

00 

2:9561636 

0.915E 

00 

3. 2600842 

0,996£ 

00 

3:164u820 

0.1 07E 

01 

2:7646163 

O.li^E 

01 

2:3995907 “ 

0.1T3E 

01 

?:3777688 

0.13CF 

01 

2'.9i6i577 

Of l38E 

01 

3. 6750324 

0.146E 

01 

2.6697096 

0.1536 

01 

1.4245761 

0.161E 

01 

0.94212050 

0.169E 

01 

0.85623063 

0.176E 

01 

1.0356017 

o.i8<e 

01 

1.5603776 

0.191E 

01 

1:2546518 

0.199E 

01 

0.88922694 

0.207E 

01 

0*V5990515 ■ 

0.2KE 

01 

0.82202307 

0,2??E 
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figure 4.1.20B 

AUTOREGRESSIVE MODEL 


MAP- 205 PASS-97 MODE- I DAS-2 
(FILTERED) 
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FICURE 4.1.21 

MAP-20S PASS 97 MODE 1 DAS-3 
ALTIMETER RANGE RESIDUALS 
POINTING * .75* 
(FILTERED) 





FIGURE 4. 1.22 

MAP-205 PASS 97 MODE 1 OAS-3 
POWER SPECTRAL DENSITY 
POINTING * .75* 
(FILTERED) 




FIGURE /k.1.24 

AUTOREGRESSIVE MODEL 


MAP-205 PASS-97 HQOE-1 OAS-3 
(FILTERED) 
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FIGURE 4.1.25 
MAP-6 PASS 7 MODE 1 DAS-1 
- ALTIMETER RANGE RESIDUALS 
POINTING =1.2® 


(UNFILTERED) 








FIGURE 4.1.27 
MAP-6 PASS 7 MODE 1 DAS- I 
AUTOCORRELATION FUNCTION 
POINTING = 1.2* 
(UNFILTERED) 
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FIGURE 4.1.29A 
MAP-6 PASS 7 MDOE 1 DAS-2 
ALTIMETER RANGE RESIDUALS 
POINTING = 1.2* 


(UNFILTERED) 




FIGURE 4.1.29B 
MAP-6 PASS-7 MODE-1 DAS-2 
altimeter range residuals 

POINTING = .2* 

CFTI TFRFO) 















FIGURE 4.1.30A 
MAF>-6 PASS 7 MODE 1 DAS-2 
POWER SPECTRAL DENSITY 
POINTING =1.2* 
(UNFILTERED) 
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FIGURE 4.1.31A 
MAP-6 PASS 7 MODE 1 OAS-2 
AUTOCORRELATICIN FUNCTION 
POINTING =1.2* 
(UM^ILTERED) 
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PIGURE 4.1.32A 

AUTOREGRESSIVE MODEL 
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FIGURE 4.1.32B 

AUTOREGRESSIVE MOCCL 

MAP-6 PASS-7 MOOE-1 DAS-2 
(FILTERED) 
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FIOURE 4.1.36 
MAP-6 PASS 7 MODE 1 DAS-3 
AUTOCORRELATION FUNCTION 
POINTING « 1.2* 
(U»riLTERED) 
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4.2 Hode 3 R— ults 





From Table 2.1 the three Mode 3 data acquisition aubmodes are 16.64 seconds* 
26.0 seconds and 101.92 seconds long. Subnode DAS-1 is a shortened replication 
Of Mode 1 OAS-l. During DAS-2 only one of the transmitted pulse pairs is used 
for deriving range information making the measurement process essentially the 
same as In DAS-l. DAS- 3 features a narrower pulse (20 nanoseconds) than DAS-1 
and DAS-2 (100 nanoseconds). 

Tabulation of the Mode 3 statistics (Table 4.3) shows little change in the 
noise sigmas for successive subnodes at low* moderate and high pointing angles. 
Examination of the power spectral plots for low and moderate pointing angles 
(Figures 4.4.2* 4.4.6* 4.4.10* 4.4.14* 4.4.18 and 4.4.22) shows that the power 
spectra are relatively flat %d.th occasional peaks which are several Db above 
the aabient level. The high pointing angle (Map 12) power spectra (Figures 4.2.26* 
4.2.30 and 4.2.34) do not appear to have as much attenuation of the higher fre- 
quencies as do the high pointing angle spectra from Mode 1 and Ifode 3. 
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TABLE 4.3 

MODE 3 DETAIL ANALYSIS 
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FIOAE 4.2.3 

MAP-45 PASS-24 MODE-3 DAS-1 
AUTOCORPELATICM FUNCTION 
POINTING s .6* 


(UISFILTERED) 










FIGURE 4.2.4 

AUTOREGRESSIVE MODEL 
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0|306E 01 
0.314E 01 
0.322E 01 
0|329E 01 
0.337E 01 
0|349E 01 
0i392E 01 
0,360E Cl 
0|368E 01 
0,375«- 01 
0(303E 01 


0,3269927 - 

0,76874942 — 
0,26366964 - 
0,17601904 
0,1794998? -- 
0,24093829 - 

0,43612812 
0,91891174 — 
1,6193213 - 
2,0623191 — 
2,2926787 - 

1,9747997 - 
1,9492809 - 

1,4356933 

1,7103440 - 

1,9937348- 
1,9401869 - 
1,1084491 -- 
1,1004876 - 

-1,7071613 

2,8872611 
■ 1,9284210 — 
0,72289829 - 
0,94801126 
0,68436982 -- 

- 1,9623203 — 

- 3,4999357 - 

l,l^^4961- 

0,51897096 

0,41870036 — 
0,93069404 — 
0,96234318 — 
1,3296819 — 
0,83633698 
0,60097696 
0,68117098 
1 

3,8246379 
1,6841681 
0,79093184 -- 
0,99906898 
0,78771436 
1,9097300 
1,9496876 
0,71833216 

0,44894629 

0,42797090 ' 

0,62784732 
1,9169923 
3,4931861 
1,7469986 







FIGURE 4.Z.T 
PASS 24 MCDC 3 D4S-2 
AUTOCCRRELATICM FUNCTION 
POINTING « .«* 
(FILTBIED) 









FIGURE 4.2.8 

AimWEGRESSIVE MODEL 

MAP-48 PASS-24 MOOE-3 OAS-2 
(FILTERED) 


GAIN 

FREQUENCY (FRACTIONAL) 


PLOT OF 

POWER TRANSFER FUNCTION 08 


0 4 8 12 


0# 

— 

0 1 12494444 


0t7C6E< 

•01“ 

— OV 14018290 


0,153E 

00 

0i28941294 

... 

0|230E 

00 

0i7162797B 


0,306E 

00 

2,8749257 

... 

0|36nE 

00 

2,1020014 


0,4aOE 

0,936E 

00 

1,3973060 


00 

^ 4972074 


C«A13E 

00 

1,8695114 

.. .. . , 

0.669E 

00 

1,7962593 



0,766E 

00 

1,5312615 

— - 

0,M3E 

00 

- 1,9076977 

— 

0,919E 

0,9«6E 

00 

00 

3,4947104 


0.2247021' 

- 1,5742883 

1,1708628 



0,107E 

0,I19E 

01 

01 


0.123E 

01 

1,6095150 


e,l30E 

01 

4,2337775 


OtlSBE 

01 

3.7255241 


0.146E 

0.193E 

01 

01 

1.2905148 

0*80926744 


0il6iE 

01 

0*91578933 

’ 
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01 

' 1,2994570 


0«176E 

01 
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* " ■ 
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01 

0*37089466 


0,191E 

01 

- 0*23454093' 


e,l9»E 

01 

0*21995119 

— 

0,(07E 

01 

0*30024219 

. . . 

0,21«E 

01 

0*7202)415 

... . 

0,222E 

01 

3,2216980 

-• *- ■ 

0|230E 

0,237E 

01 

3,0598107 
— " 1,8457403’ 


01' 


P,2<5£ 

01 

1,1907998 

..... 

e.293E 

01 

1,1905180 

- 

0,26gE 

01 

0,87761220 


0,2*>t 

01 

0.57760281 


C.?' 


0.49617847 


l,283E 

01 

0,60671124 


0,291E 

01 

0,92296423 


8.299E 

01 

- 0,903027'»0 

.... 

0,39BE 

01 

0.99657035 


0,914E 

01 

0*42358647 


0,322E 

01 

0*40866757 

• 

0.329E 

01 

0'*94»‘86612 


e,337E 

01 

2,0665225 

* 

0«349E 

01 

1,6620114 

.. 

0.392E 

01 

0,66903151 


8.S60E 

01 

0,46094043 


B,3f»E 

01 

0,40174231 


0.379E 

01 

0,702*4807 


0,3B3E 

01 

1,3397390 






4-57 



FIGURE 4.2.9 

MAP-45 PASS 24 MODE 3 OAS-3 
altimeter range RESIDUALS 
POINTING 3 .6* 

(FILTERED) 








FIGURE 4.2.10 

MAP-45 PASS 24 MODE 3 OAS-3 
POMER SPECTRAL DENSITY 
POINTING » .6* 
(FILTERED) 



















FIOJRE 4.2.12 

ALfTQREGRESSIVE HDOEL 


MAP-4S PASS-24 MOOS-3 OAS-3 
CFILTBCO) 


GAIN 


PLOT OF 


FREQUENCY 

(FRACTIONAL) 

0. 

0.766E< 

'Oi 

»,.*97359UE 

0.82825438E 

0.193E 

00 

C.14!»17544 

0.231E 

00 

0.4'»«222!'9 

0.306E 

00 
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0.38TE 

Ou 

1.1039$46 

0.469E 

00 

0.^1053249 

0.536E 

00 

U. 50493929 

0.413E 

Ou 

0.9J434K89 

0.449E 

Ou 

2.4443677 

0.766E 

00 

2,4744278 

0.843E 

CO 

1.3tUM3 

0.919E 

Ou 

1.0144^^1 

0.996E 

ou 

1,0324759 

0.197E 

01 

fc.97783Hll 

O.U!>E 

01 

0,79393343 

0.12SE 

01 

C.72052C26 

0.13UE 

01 

C.07V32554 

0.13eE 

01 

1.476?1?5 

0.14CE 

01 

2,2305469 

O.l^SE 

01 

1,4047605 

0.161E 

01 

1,1039624 

0.169E 

01 

1, 0927495 

0.176E 

01 

1,5260647 

o.ie^E 

n 

2.0444108 

0.J91E 

01 

1,5344932 

0.199E 

01 

1.0170C01 

0.2C7E 

01 

0,50.»C2«'1V 

0.214E 

01 

1,1237979 

0.222E 

01 

1.74f,4435 

0.230E 

01 

2, 2124757 

0.237E 

01 

1,4627775 

0.245E 

01 

1,2421227 

0.233E 

01 

1.C«'30372 

0.26.IE 

31 

1, 11225978 

O.E^oE 

01 

0.93264609 

0.276E 

01 

11.0314341 4 

0.283s 

01 

0.79007713 

0.?9lE 

01 

0.62401C3Q 

0.299E 

01 

t. 69438570 

0.3C4E 

01 

0.52529007 

0.3l«E 

01 

f .91017372 

P.3??E 

Cl 

C.5H494#42 

0.374F 

05 

i,'»162.320 

0.33’l 

01 

1.1164;?5 

0.9<>JS 

C} 

l,irp.R744 

0.9*/7<i 

01 

f..VV‘'VCV?7 

O.SfrCE 

01 

(i.52.f.45;»ia 

0.348E 

01 

0.(494462 

0.374E 

Cl 

1,11505005 

0.38J{= 

Cl 
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FIGURE 4.2.13 

^ MAP-10 PASS 8 P«30E 3 OAS-1 
. ALTIMETER RANGE RESIDUALS 
POINTING « .75* 
(UNFILTEREO) 








flUT0C0Ri^eLATI9M 

•0.1 0.0 e.t 




FIGURE 4.2.16 

AUTOREGRESSIVE MODEL 

MAP-10 PASS-8 MQOE-3 OAS-1 
(UNFILTERED) 


GAIN 


FREQUENCY 

(FRACTIONAL) 

0, 


1,0105368 

d,766E< 

•01 

1,1067072 

0,1!>3E 

80 

1,4641976 

0.230E 

00 

2,1825912 

0,30fE 

00 

2.4235699 

0.3B3E 

00 

1,6625982 

0.^63E 

00 

1,2250007 

0.536E 

?c 

1, 30*3572 

0.613E 

00 

2.17(i7te4 

0.669E 

00 

2.8862337 

0,766E 

00 

1.1727360 • 

0,84?E 

00 

0,54713986 

0,9196 

00 

( .39446134 

0.99AE 

00 

(..•.2455505 

0.1076 

01 

C. 69941184 

0.1156 

91 

1 ,6692671 

0.123E 

01 

2,5563606 

0.l3uE 

01 

J. 5209658 

O.lSBE 

01 

1,2169028 

tf,l<6E 

01 

J ,5138198 

0.153E 

01 

1.9298571 

0.1616 

01 

1 ,3161300 

0,1696 

01 

c.03J68s:i7 ■ 

0.1766 

01 

5, 73763500 

0.104E 

01 

(,.8 7355091 

0.191E 

01 

C.(8303C16 

0.199E 

01 

C. 55023374 

0.2076 

91 

0.34720650 

0.2146 

01 

C. 29/59674 

0.2226 

01 

f..37c.70M9 

D.230E 

01 

C.7C)79t98 

0,2376 

01 

1.56S2C76 

0,2<5E 

01 

1.2^253)0 

0.253E 

01 

( .74241608 

0.260E 

91 

t. 73065197 

0.26PE 

01 

1,2490515 

0.2766 

01 

2.3U9?.?? 

0.2S36 

91 

1,1177405 

0.291E 

91 

6.53877594 

0.299E 

91 

0.42H57448 

C.3C66 

91 

0.54725443 

0.3146 

01 

1,1109218 

0.3226 

01 

1,9953428 • 

0.329E 

91 

0.93956354 

0.337E 

91 

l).5lV05r.64 

0.3456 

01 

0.45542252 

0,3526 

91 

(•.6l5l0l'J6 

0.36CE 

01 

1,3254243 

C.36eE 

01 

3,1321541 - 

0.375E 

01 

l,7*)?25e6 

0,383E 

91 

r.898C4231 


PLOT OF 

POMER TRANSFER FUNCTION DB 



FIGURE 4.2.17 
MAP-10 PASS 8 MODE 3 OAS-2 

ALTIMETER RANGE RESICXMLS 

POINTING » .75* 
(umLTERED) 




p-30.0 
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FIGURE 4.2.19 
MAP-IO PASS 8 MODE 3 OAS-2 
PONER SPECTRAL OENSITV 
POINTING = .75* 
(UTFILTERED) 












FIGURE 4.2.20 

AUTCREGRESSIVE MODEL 

MAP-10 PASS-8 MOOE-3 DAS-2 
(UNFILTEREO) 

GAIN 


frequency 

(FRACTIONAL) 

0. 

0|TA6E« 

’01 

2,1694376 

ll9635978 

0.193E 

00 

l‘.$087603 

0i230E 

no 

1',099.<19i 

o.soen 

00 

Ct08632Sl2 

0.3A35 

00 

0.86857620 

0.46pE 

00 

i;07ci932 

0.536E 

00 

1.5556476 

0.613E 

00 

1. 914414? 

0,489E 

00 

1'.497320i 

0.766E 

00 

1.07449C6 

0,A43F 

CO 

0.94043324 

0.919E 

00 

lln5327Jo 

0.996E 

00 

1'.4377652 

0.107E 

01 

1'.9o2<607 

0.115E 

01 

1;6!>86573 ' 

0.123E 

01 

1^0033259 

0.130F 

01 

0.75803253 

0.J38E 

01 

0.64657454 

0.146E 

01 

0. 70381793 

0.193E 

01 

l'. 0394666 

0.161E 

01 

2;2927987 

0.169E 

01 

5.0393634 

0.176E 

01 

2.4656031 

0U84E 

01 

1'.1555162 

0.191E 

01 

0.62806005 

0.199E 

01 

0.02955968 

0.207E 

01 

1.0224749 

0.214E 

01 

1.2108283 

0.??2E 

01 

l'.l02l98l 

0.2S0E 

01 

0.900l90(*S 

0.237E 

01 

0.80742226 

0,?45E 

01 

0.80364373 

0,2$3E 

01 

0.80796835 

0i260E 

01 

0.77667053 

0.268E 

01 

0.765948(6 

0.276F 

01 

0.8)285827 

0.283E 

01 

0.84530048 

0.291E 

01 

0.71054571 

0.299E 

01 

0.51186176 

0.304F 

01 

0.39799834 

0.314E 

01 

0.3764966»,' 

0.322E 

01 

0.44132113 

0.329F 

01 

0.59893347 

0.337r 

01 

0.76249796 

0.345E 

01 

0.75856643 

0.3^2H 

01 

0.69639123 

0,34nE 

01 

0. 72828519 

0 « 3 AdF 

01 

0.68957006 

0.37iF 

01 

1.0959658 

0.393E 

01 

1;1384738 
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8 
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FIGURE 4.2.21 

MAP-10 PASS S MCXIE 3 OAS-3 

ALTIMETER RANGE RESIDUALS 

POINTING * .75* 

IFlLT Crg )) 






CCO*?RELfiriON 



FIGURE 4.2.23 
MAP- 10 PASS 8 MODE 3 DAS-3 
AUTOCORRELATION FUNCTION 
POINTING » .75* 
(FILTERED) 


FIGURE 4.2.24 

AUTOREGRESSIVE MODEL 

MAP-10 PASS-8 HOOE-3 OAS-3 
(FILTERED) 


FREQUENCY 

(FRACTIONAL) 

0. 


0.17293816 

Of 766E" 

■01 

0.20897720 

0.153E 

00 

0.39119275 

Oo23uE 

00 

ll424ti3o 

C.306E 

00 

310996371 

0.393E 

00 

0.t)3363l7o 

0.460E 

00 

0.45366816 

0.536E 

00 

0.41757502 

0.613E 

00 

0.5768o5?0 

0.689E 

00 

l’. 1393165 

0.766E 

00 

2‘.268?577 

0,843g 

00 

210387822 

0.9)9E 

00 

Il402l38l 

0,996E 

00 

ll?37686i 

0.107E 

01 

113387306 

0.115E 

01 

115051104 

0.123E 

01 

115633835 

0.13CE 

01 

1.5708223 

0.133E 

01 

ll6436l93 

0.146E 

01 

116925624 

0.153E 

01 

1.4879822 

0.161E 

01 

111206597 

0.169F 

01 

0.84047654 

0.176E 

01 

0.69886602 

0.184F 

01 

0.67325359 

0.l9lE 

01 

0, 76941682 

0.199E 

01 

ll0640303 

0.207E 

01 

i; 7363494 

0.214E 

01 

2.5511/75 

0.222E 

01 

211707757 

0.230E 

01 

114597956 

0.237E 

01 

ll 1424450 

0.245.F 

01 

ll 1025326 

0.283E 

01 

ll2l87463 

0.260E 

01 

1 12867927 

0.263E 

01 

110807685 

0.276E 

01 

0.77127442 

0.2B3E 

01 

0.55633029 

0.291E 

01 

0.4<133502 

0.2?9E 

01 

0,39479892 

0.3C6F 

01 

0.40224777 

0.314E 

01 

0.46877608 

0,372E 

01 

0.61451277 
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0.329E 01 
0.337E 01 
0.315F 01 
0.S52E 01 
0.360F 01 
0.36BT 01 
0.375E 01 
0.383E 01 


0.b39!>(i616 
i;03l60'*8 
ll073967p 
1;0777434 
1.1708232 
1‘.3S4393S 
1 *. 4668861 
1.3927333 
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FIGURE 4.2.26 


MAP-12 PASS-9 MDOE-3 DAS-1 
POWER SPECTRAL DENSITY 
POINTING = .6* 
(UNFILTERED) 




FIGURE 4.2.27 
MAP-12 PASS-9 MODE-3 DAS-1 
AUTOCORRELATION FUNCTION 


POINTING = .6 
(UNFILTEREO) 





FIGURE 4.2.28 

AUTOREGRESSIVE MODEL 

MAP-12 PASS-9 MOOE-3 OAS-l 
(UNFILTERED) 


FREQUENCY 

(HZ) 


0 . 

Ot 766 Es01 
0.153E 00 
o.23ge oa 
0.306E 00 
0t3R3E 00 
0.460E 00 
0.536E 00 
0.^13E 00 
0.689C 00 
0.766E 00 
0.843E 00 
Oi’l’E 00 
0.996E 00 
0*J07E 01 
0ill5E 01 
0>123E 01 
0il30£ 01 
0.138E 01 
0.146E 01 
0 i 153E 01 
C.161E 01 
0«169E 01 
0il76E 01 
0.184E 01 
0.191E 01 
0,199E 01 
0.207E 01 
0»2l4E 01 
0.222E 01 
0;Z30E 01 
0.237E 01 
0.2«!>E 01 
0.253E 01 
• 0.260E 01 
0.268E 01 
0.276E 01 
0|283E 01 
0;29lE 01 
0.29VE 01 
0«308E 01 
OiSlAE 01 
0|322E 01 
0'i 329E 01 
0i33/E 01 
0.345E 01 
0.382E 01 
0i3A0e 01 
0.368E 01 
9(37i>E 01 
0»363E 01 


GAIN 

(FRACTIONAL) 
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3.121H798 
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1.2305630 
1.7949303 
3.66686<2 
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1.9296267 
1.4495126 
1.5340729 
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1.1154949 
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0.87169117 
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0. 40116737 
0,39363704 
0. 45o35459 
U. 51710384 
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FIGURE 4.2.29 

MAP-12 PASS 9 MGOe 3 DAS-2 
- ALTIMETER RANGE RESIDUALS 
POINTING a l.l* 
(UNFILTEREO) 
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FIGURE 4.2.31 

MAP-12 PASS 9 MCX3E 3 OAS-2 
AUTOCORRELATION FUNCTION 
POINTING =1.1* 


(UNFILTERED) 










FIGURE 4.2.32 

AUTOREGRESSIVE MODEL 

MAP- 12 PASS-9 MOOE-3 OAS-2 
(UNFILTERED) 


GAIN 

FREQUENCY (FRACTIONAL) 
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POWER TRANSFB? FUNCTION OB 
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31 
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01 
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31 

0.?(t-E 
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91 
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F1€UNE 4.2.33 
MAP-12 PASS 9 MODE 3 DAS-3 
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FIGURE 4.2.36B 

AUTOREGRESSIVE MODEL 


MAP- 12 PASS-9 MOOe-3 DAS-3 
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4.3 Mode-5 Results 


From Table 2.1 the three Mode 5 data acquisition submodes (DAS-1, MS-2 and 
MS-3) are 16.64 seconds, 102.96 seconds and 53.04 seconds long. DAS-1 is a 
shortened replication of Mode 1 DAS-1. MS-2 uses pulse compression techniques 
to increase range measurement accuracy. DAS-3 transmits a 20 nanosecond pulse 
for comparison with the DAS-2 pulse compression results. Discussions with NASA/ 
Wallops have Indicated that, out of the three Mode 5 data sets (Maps 47, 74 and 66) 
selected for detailed analysis the pulse compression was operating only during 
Map 74. 

Examination of Table 4.4 shows that for low, moderate and high pointing 
angles (Maps 47, 74 and 66, respectively) the noise sigma is largest for sub- 
mode DAS-1 followed by DAS-2 and DAS-3. Table 4.5 shows the ratios of the 
DAS-1 and DAS-2 sigmas to the DAS-3 sigma. Note that for DAS-1 this ratio in- 
creases from 1.58 for low pointing angles to 2.67 for high pointing angles 
(Map 66) while for DAS-2 this ratio is constant at about 1.3 for all pointing 
angles. 

At low pointing angles (Map 47) the DAS-1 power spectrum (Figure 4.3.2) and 
autocorrelation function (Figure 4.3.3) do not indicate any significant deviation 
from a white noise process. The DAS-2 power spectrum (Figure 4.3.6) shows con- 
siderable attenuation of frequencies higher than 1.25 Hz. The autocorrelation 
function (Figure 4.3.7) shows 8 of the first 20 autocorrelation coefficients 
exceed the plus or minus two sigma limits for a white noise process. The DAS-3 
power spectrum (Figure 4.3.10) also shows some attenuation of the higher fre- 
quencies and 5 of the first 20 autocorrelation coefficients exceed the two sigma 
limits for a white noise process. 

At moderate pointing angles (Map 74) the DAS-1 power spectrum (Figure 4.3.14) 
is relatively flat and only two of the first twenty autocorrelation coefficients 
(Figure 4.3.15) exceed the two sigma limits for a white noise process. The DAS-2 
power spectrum (Figure 4.3.18B) sho;Ts 10 to 15 db attenuation at frequencies 
higher than 1.25 Hz and the DAS-3 I'owe." spectrum (Figure 4.3.22) shows a similar 
attenuation at frequencies higher than 1.0 Hz. Comparison of the DAS-2 and 
DAS-3 autocorrelation functions (Figures 4.3.19 and 4.3.23) show that the 
MS-3 has a higher degree of serial correlation. 


TABLE 4.5 
MODE 5 

RATIO OF SUBMODE 

NOISE SIGMAS TO OAS-3 NOISE SIGMA 

SIGMA/ 


MAP 
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DAS- 3 SIGMA 

47 

DAS-1 

1.58 



1.36 

74 
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1.93 


( 1 )DAS-2 

1.30 

66 
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2.67 


(2)DAS-2 

1.30 
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The effect of high pass filtering on the Map 74 DAS-2 data is shown by 
comparing the unfiltered and filtered time series plots (Figures 4.3.17A and 
4.3.17B) and power spectra (Figures 4.3.18A and 4.3.18B). The unfiltered time 
series plot shows a well defined slope of about .2 meters per second for the 
first 35 seconds of data. This produces a large low frequency (0 to .10 Hz) 
spike in the corresponding power spectrum (Figure 4.3.18A). Conq>arlson of the 
unflltered data power spectrum (Figure 4.3.18A) and the filtered data power 
spectrum (Figure 4.3.18B) shows that filtering had no appreciable effect on the 
spectral coefficients for frequencies greater than .2 Hz. 

At high pointing angles (Map 66) the time series plots (Figures 4.3.25, 

4.3.29 and 4.3.33) show that the measurement noise is domin.. ted by low frequency 
( .2 Hz) oscillations. The power spectra (Figures 4.3.26, 4.3.30, 4.3.34) show 
considerable attenuation of spectral coefficients at frequencies greater than 
.25 to .5 Hz. Note that the DAS-2 and DAS-3 spectra are plotted at 10 Db per 
division instead of the 5 Db per division used for all other power spectral plots. 
The autocorrelation function (Figures 4.3.27, 4.3.31 and 4.3.35) are smooth 
damped oscillation type curves. 
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FIGURE 4.3.1 

MAP-47 PASS 25 MODE 5 DAS-1 
ALTIMETER RANGE RESIDUALS 
POINTING = .5* 
<Ut>FILTERED) 























FIGURE 4.3.3 

MAP-47 PASS 25 MODE 5 DAS-l 
AUTOCORRELATION FUNCTION 
POINTING * .5* 
(UNFILTERED) 
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FIGURE 4*3.7 


MAP-47 PASS-25 MQDE-5 DAS-2 
AUTOCORRELATION FUNCTION 
POINTING = .5® 
(FILTERED) 



FIGURE 4.3.8 
AUTOREGRESSIVE MODEL 
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FIGURE 4.3.12 

AUTOREGRESSIVE MODEL 


MAP-47 PASS-5 MOOE-5 DAS-3 


(FILTERED) 


frequency 


GAIN 

(FRACTIONAL) 


PLOT OF 

POWER TRANSFER FUNCTION DB 


0 . - 

-3. 766E«i)l 

- Q,153F 00 
0,23Cfe 00 

--9v306s -00- 
-0,3636 00- 
-0,'«606 )u- 
0.t>36fc OO 
■ 0.fl3E 00 

- 0.6806 09- 
— 0 .^ 666 — 99 — 

- 09 - 
--0.919= 09- 
-0,9966 00 

- 0;1076 01 - 

01 - 

■0.-1236 01 - 
0,1306 01- 
0,1 35E 01- 
0.146= 01 - 
0;1536 01 

0,1616 Cl 
0.1 696-01- 
0.176c 01 

0,194fc 01 
C;191E 01 
0.199E 01 
0.207k 01 
0.214E- 01— 
0,2226 01 
C.2306 Cl 
0.237E 01 
0.245": 01 
0,2536 01 
0.2606- Ol" 
0.2636 31 
0.2766 01 
0,2536 01 
0,2916 01 
0.2996 01 
0,5066 01 
0j3J.<6 01 
-0.-3226-01- 
0,3296 01 
0,3376 01 
0.340c Oa 
0.3526 01 
0.3!';Gc 01 
9, 3<> = E. Oi- 
0,375c 01 
0 . 343 = 01 


-- C.V6439741F-01 
0.11197408 
C.l5tf2i.845 
C. 44870929 

2.325I.V17 

“ 6,3226374 

— 2.220)411 
--1.93l"316 • 

- 3,023»703 

-■4,047‘'5V3 - 

-2,6511321- 

--2,C92«:234 ... 

-- 2,1234053 

- 3,5463643 

■ 2, 0831431 

— r-.'Si>9t«<ofi3 

- (i.9c5n365 

- 1,5275916 

- I,rl7)587 

- C.9t2rl276 
— CTT‘t4 20f)61 

--1,179702? 
1.79C6929 
1, 55^0349 
• C. 95811410 
C.74«C«479 
— C77219291? 

- I,l955l9§ 

■ l.o-l2o33 
1. 3227587 

t. 9^759750 
C. 93205758 

- -i;.t-»0446i - 

1,.*9';3960 
1,5363762 
1,0152435 
(1.69291*85? 

1‘. 57876002 

f. .579sn;;o0 
0,6-»845jHl 
-C-.-727014i.5- 
C . 73lt032«»6 
t , •523-(5ft i4 
0,41>2.-,6275 
C.5763(WJ72 
C. 322 334 96 

- V. >1802201 - 

C.3/12‘»‘il4 
0,<945O;j9 


OBIGIHAL 
OP POOR QUALfPP 






FIGURE 4.3.13 

MAP-74 PASS 39 MODE 5 DAS-1 
ALTIMETER RANGE RESIOUALS 
POINTING = .85“ 
(UNFILTERED) 







FIGURE 4.3.14 

MAP-74 PASS 39 MODE 5 DAS-1 
POWER SPECTRAL DENSITY 
POINTING = .85* 


ItFFILTERED) 




FIGURE 4.3.15 

MAP 74 PASS-39 MOOE-S DAS 1 
AUTOCORRELATION FUNCTION 
POINTING s .85* 


(UNFILTERED) 
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FIGURE 4.3.176 
HAP-74 PASS 39 MOOE 5 DAS-2 
ALTIMETER RANGE RESIDUALS 
POINTING = .85* 
(FILTE R ED) 








FIGtjRE 4.3.1S8 
MAP-74 PASS-39 MOOE-5 OAS-3 
POWER SPECTRAL DENSITY 
POINTING = .05* 


FILTERED) 











FIGURE 4.3.20 

AUTOREGRESSIVE MODEL 

MAP-74 PASS-39 MQOE-S OAS-2 
(FILTERED) 
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FIGURE 4.3.21 

MAP-74 PASS 39 MODE 5 OAS-3 
ALTIMETER RANGE RESIDUALS 
POINTING s .8S* 
IFILTERED) 




FIGURE 4.3.22 

MAP- 74 PASS 39 MODE 5 OAS-3 
POWER SPECTRAL DENSITY 
POINTING s .85* 
(FILTERED) 











FIGURE 4.3.23 

MAP-74 PASS 39 MODE 5 DAS-3 
AUTOCORRELATION FUNCTION 
POINTING = .85* 
(FILTERED) 








FIGURE 4.3.24 

Al/rOREGRESSIVE MODEL 


MAP- 74 PASS-39 MODE-5 DAS-3 




(FILTERED) 
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FIGURE 4.3.25 
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FIGURE 4.3.27 
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As discussed In Section 2.2, the coefficients for an autoregressive model 
of order 20 are computed from the autocorrelation function. This model is use- 
ful in that it describes the serial correlation structure of the time series 
with a small set of parameters. The analytical power transfer function computed 
from the autoregressive model is actually a smoothed representation of the power 
spectral density of the original time series. 

For Map 39 DAS-1 comparison of the power spectrum (Figure 4.1.2) and the 
autoregressive model power transfer function (Figure 4.1.4) shows that both 
curves have the same overall shape although the power spectrum, which has 108 
coefficients, shows more detail than the power transfer function, which was 
derived from 20 autoregressive coefficients. Table 4.4.1 compares the frequencies 
of major power spectrum peaks with peak frequencies of the power transfer func- 
tion. Peaks in the power spectrim separated by less than .15 Hz were merged into 
one peak in the po~rer transfer function. Power spectrum peaks separated by .3 
Hz are well defined in the power transfer function. 

For Map 45 DAS-3 Table 4.4.2 shows a similar comparison of peak frequencies 
for the power spectrum (Figure 4.2.10) and the autoregressive model power trans- 
fer function. Again, the power spectrum, which in this case has 215 coefficients, 
shows much greater detail than the power transfer function which is derived from 
20 parameters. Peiiks separated by less than .12 Hz were combined into a single 
peak in the power transfer function. 

Overall qualitative evaluation indicates that with 20 autoregressive coef- 
ficients the autoregressive model power transfer function provides a smoothed 
representation of overall power spectrum shape and will reproduce individual 
peaks which arc separated hy at least .3 Hz. 




TABLE 4.2.2 
MAP 45 DAS-3 

LOCATION OF POWER SPECTRUM PEAKS 


POWER SPECTRUM AUTOREGRESSIVE MODEL 

(FIGURE 4.2.10) TRANSFER FUNCTION 

(FIGURE 4.2.12) 
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1.0 INTRODUCTION 


The Residual Analysis Program (RAP) is designed to perform a si. . Istical 
analysis of altimeter residual data. The program performs data screening, 
displays various descriptive statistics and plots, tests for stationarity and 
normality and computes the power spectral density and autocorrellogram. Each 
input data set is assumed to contain 3 colvmms of data which are: 

Column 1 - Altimeter Resldtial data 

Column 2 - Other time series data sequence 

Column 3 ^ Time labels corresponding to entries in Columns 1, 2. 

The specific operations performed by RAP on Columns 1 and 2 of the data 
array include the following: 

1. Computation of overall mean, standard deviation, maximum, minimum, 
range and "t" statistic; 

2. Time series and histogram plots; 

3. Partitioning of each column into subgroups and computation of subgroup 
means and standard deviations; 

A. Computation of statistics to assess the homogeneity of subgroup 
variances, means and the underlying probability distributions; 

5. Computation of statistics to test the hypothesis that the residuals 
follow a normal probability distribution: 

6. For each column, computation of the power spectral density and 
autocorrelation function and autoregressive model coefficients. 

7. Computation of the cross correlation function and cross spectral 
density for columns 1 and 2 of the data array. 

RAP contains logic for the ex^ractlon of data from an input tape between 
given start and stop dates and times. Outlying observations in either Column 1 
or Column 2 of the data set can be deleted if the absolute value exceeds a 
user specified edit threshold. Deleted data points are replaced with the error 
code - 7777. and are ignored except in the computations of power spectral density, 
autocorrelation coefficients, cross correlation coefficients and cross spectral 
density. For these computations the deleted data points (error code - 7777.) 
are replaced with a local average. 


A-3 


2.0 DESCRIPTIVE STATISTICS 


RAP outputs various descriptive statistics and plots designed to give the 
user a feel for the underlying data structure. The sample size, for each 
data column is defined as the total number of data points minus the number of 
edited data points. 

2.1 Time Series Plot 

The time series plot is simply a display of the raw input data (after edit- 
ing) in time ot ier. The abscissa (x axis) contains the data time tags a.id the 
ordinate (y axis) is the value of the associated data point. 

2.2 Histogram 

The histogram is a means of graphically representing the sample probability 
density function. The abscissa is divided into equally spaced intervals and 
the ordinate is scaled to give the percentage of data values falling within each 
Interval. 

2.3 Sample Mean 

The sample mean x, of a variable x, is the average value of the N 
measurements of that variable. This statistic is a measure of the centrality 
or location of the distribution. It is calculated from the expression 


X 


N- 

I X. 

i-0 


( 2 . 1 ) 


2.4 Sample Standard Deviation 

The sample standard deviation, s, is a measure of the dispersion of the 
N variable measurements about the sample mean value x . As the mean x 
corresponds to the center of gravity of the mass of the distribution, M times 
the moment of second order (v-iriance) corresponds to the moment of Inertia of 
the mass with respect to a perpendicular axis through the center of gravity. 

The standard deviation, s, corresponds to the positive square root of this 


variance so as to keep the statistic In the units of the variable. 

, , N-1 , 

“ ■ jo 

2 

The reason for dividing by (N-1) Is so that s Is an unbiased estimator of 
the population variance. 


2.5 Standardized Variates 

Prior to computing the power spectrum, autocorrelation function, cross corre- 
lation function, cross spectrum and autoregressive coefficients described in Sec- 
tion 2.6-2.10, all edited data points are replaced by a local average defined as 

®i“ 2^^i'^^i-l^ 


where f^^ and are those two non-edlted data points directly preceding the 

edited data point. This replacement Is necessary because the algorithms used 
In computing the power spectrum and autocorrellogram assume that the data points 
are all piesent and equi-spaced in time. After this replacement the mean and 
standard deviation are computed as 


X • 


1 

N 


N-1 


I 

1-0 


(2.4) 



N-1 

(x . - x) 

1-0 


(2.5) 


Note that for these formulas N Is the total number of time points in the data 
set since all edited points have been replaced with a local average. The mean, 
Eqn. (2.4), is exactly as defined In Eqn. (2.1) but the standard deviation, 

Eqn. (2.5), contains a divisor of N Instead of N-1 as in formula (2.2). 

This is done so that the computed autocorrelation coefficients will be as defined 
by Blackman and Tukey (Reference 1, pp. 53). 
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The sti^ndardlzed time series variates are computed as 







( 2 . 6 ) 


In the following descriptions of power spectrum and autocorrelation computations 
it Is assumed that the edited data points have been replaced and the resultant 
time series standardized as per Eqn. (2.6). 


2.6 Power Spectral Density 

For the discrete time series (Xj, j«0,l,2, • • *N-1) of N equl-spaced data 
points the Fourier transform is defined as 


N-1 


(2.7) 


where 1 ■ / -1 and the {a^^} are In general complex valued. 

If At Is defined as the time Interval between data points then the total 
interval Is given by 

T - NAt (2.8) 


The frequency resolution Is Inversely related to the total time interval as 


Af ■ 


1 . 

T NAt 


(2.9) 


Thus the frequency associated with the Fourier coefficient is given by 


f^-Uf 


( 2 . 10 ) 


If the time series (Xj) 
RAP» the Fourier coefficients 


is real valued, and N is even, as assumed in 
and complex conjugates. Thus 
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if l8 expressed as 


0^ - k-0,1.2 N-1 


( 2 . 11 ) 


Then Is given by 


- a^-ib^ - \ k-l,2...N/2-l (2.12) 

it 

where signifies complex conjugate. The coefficients <“0 

unique. 

In RAP the Singleton fast Fourier tr.>nsform (FFT) algorithm* described in 
Reference 8, is used to compute the complex Fourier coefficients defined in 
Eqn. (2.11). The selection of this algorithm was based on the competitive 
evaluation published in Reference 5. 

The power spectral coefficients are defined as 


Pk • 2N(a^+b^) k.l,2---(N-l)/2 

^N/2 “ **^®N/2‘‘''’n/2^ 


(2.13) 


Eqns. (2.11) and (2.12) demonstrate that the power spectral coefficients for 
k ■ N/2+1,"*N~1 are mirror Images of the spectral coefficients for k-1,2 
N/2-1. Thus a time series containing N data points yield N/2+1 unique 
spectral coefficients. From equation (2.10) the highest frequency for which a 
unique spectral coefficient can be computed is given by the nyquist frequency 
defined as 


f„ - ( ? ) Af (2.14) 

. q 2 

This definition of the power spectral coefficients from reference 2 is program- 
med in RAP. If the time series, {x^} is a white noise process (i.e. with no 
serial correlation) then each individual power spectral coefficient 
{Pj^, k-1,2* •* Cv'-l)/2) will be distributed as a chi square variate with two 
degrees of freedom and the coefficients yill be uncorrelated with each other. 
(Reference 1) . ^ 


In general the discrete time series {x^^} will have been obtained by 
sanpllng a continuous process or by subsampllng a much denser discrete time 
series. In either case. If the original data contained components at frequencies 
higher than the power spectrum computed for {x^} will contain contribu- 
tions due to these hl^er frequencies. The contributions will be folded such 
that the power spectral coefficient of the original series at frequency 
would appear in the spectral coefficients of {X|^} at frequency f^-pdf. 

An option Is available to smooth the Fourier coefficients with a Hanning 
window prior to computing the power spectral coefficients as defined in Eqn. 
(2.13). The smoothed coefficients are given by 


• -.250j^^j^+.50j^ - k-l,2‘**N/2-l 


(2.15) 


The end coefficients and 


are smoothed as 


.5 - .5 


^/2 " “ n / 2-1 “ n /2 


(2.16) 


The power spectrum computed from the {S^} is called a modified power spectrum 
or modified perloJogram. 


2.7 Autocorrelation Coefficients 

ror a data sequence {x^, j*0,l,2> • 'N-l} the autocorrelation coefficients 
are defined as 


r 


k 


1 

N-k 


N-k-1 


I 

j-0 


XjXj+k 


(2.17) 


These coefficients can be computed directly from Eqn. (2.17) or by computing the 
power spectral coefficients of (xj) and then using an FFT algorithm to take 
the inverse Fourier transform of the power spectrum. Both methods are avellable 
in RAP. The direct method (from (2.17)) will be more efficient If only a few 
coefficients are to be computed, otherwise the second method Is preferred. 

The computational steps for obtaining all N of tne autocorrelation 

coefficients using an FFT algorithm Is as follows: 



1. Imbed the original time series of length N in a new series of length 
2N where the last N entries are zero and the first N entries are 
the original time series. 

2. Compute 2N power spectral coefficients for the new time series. 

3. Take the inverse Fourier transform of the power spectrum to obtain 
2N unsealed autocorrelation coefficients {Rj^, k»l,l* • •2N-1}. 

4. From the first N unsealed coefficients compute N scaled coefficients 
as 

Tr " \ ’ •2' ■ 

The padding with zeros in step 1. is required because of the circular Indexing 
employed within the FFT algorithm. If this padding were not utilized the resul- 
tant autocorrelation coefficients would be given by 

'k ■ 'k ^ Vk «-w) 

where r^^ are the coefficients obtained by taking the Inverse Fourier transform 
of the power spectrum of the original unpadded time series and rj^ are the 
autocorrelation coefficients defined in Eqn. (2.17). 
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2.8 Croai Spectral Coefflcltnf 

Given two real data aequencea {x^, 1>0, 1,2...N-1} and {y^, J*0» 1.2...N-1), 
the Fourier tranaforna of these data sequences are defined from Equation 2.7 as 


Y, - 


From Reference 3 the 


^ I X, exp k-0,l,2...H-l 

" i-O ^ " 

1 21l^1c 

^ I Yj exp k-0,l,2...N-l 

j-0 

cross spectral coefficients are defined as 
®k“^\ k-0,l,2...N-l 


( 2 . 20 ) 


( 2 . 21 ) 


Uhere * slgnlflas complex conjugate and X^, are in general complex. Note 
that for real data sets 

^N/2+k " ^N/2“k k*0, 1 j 2. . .N/2“l (2.22) 

In RAP the Singleton fast Fourier transform algorithm, described In 
Reference 4, Is used to compute the complex Fourier coefficients defined in 
Eqn. (2.20). An option is available, as described in Section 2.6^ to smooth 
the Fourier coefficients with a Hanning window prior to computing the cross 
spectral coefficients as in Eqn. (2.21). Because of the "mirror isiage" effect 
described in Eqn. (2.22) the cross spectral coefficients, G^, are computed 
only for k*0,l,2. ..M/2. 
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2.9 Cross Correlation Coefficients 


Given tvo real data sequences and {yj^,l«0,l,2. . .N-1}, 

the cross correlation coefficients are defined as 

||^j E *1^ ^1+k -*(N~1) » . • .0» • • «N~1 (2.23) 

where the sumDatlon over i is such that 14-k is always between 0 and N-1. 

The computational steps for obtaining all 2N'’l crosj correlation coefficients 
{Rj^} using the Singleton FFT algorithm is as follows: 

1. Imbed each time series of length N in a nev series of length 2N as follows 

{xp • {Xq,Xj^...Xj,_j^, Ojj, 0jh-1***°2N-1^ 

{yp ■ yQ» 

2. Conq>ute the 2N complex Fourier coefficients and Y^} for each of 

these padded series as in Eqn. (2.20). 

3. Conq>ute the cross spectral coefficients as 

®k “ Vk k-0,l»2...2N-l 

4. Take the Inverse Fourier transform of the power spectral coefficients {G^} 

to obtain 2N unsealed cross correlation coefficients k* 

...N-1} 


3. Compute the scaled cross correlatioii coefficients as 

N^^ k--(N-l)....0,...N-l 


(2.24) 


2.10 Autoregressive Model 


In Sections 2.5 and 2.6» the serial Interrelationships of a discrete 
stationery tlae series are described by N/2 power spectral coefficients and 
N-1 autocorrelation coefficients. Box and Jenkins (Reference 2) describe 
parametric modeling techniques which represent these serial Interrelationships 
with a much smaller nuaber of parameters. Their approach consists of modeling 
the discrete series as the output of a linear digital filter with white noise 
Input. The filter output Is given as 


where a Is the standard deviation of the white noise process Z Is the complex 
exponential £ is the frequency In Hz« T Is the sample time In seconds, 

and F(Z), H(Z) are finite polynomials in Z. The parameters to be estimated are 
the polynomial coefficients. In general this requires Iterative nonlinear 
estimation procedures. Parzan (Reference 7) and Graupe (Reference 4) suggest 
an Intermediate procedure which requires the estimation of a larger nud>er of 
parameters but requires only a single Iteration linear estimator to recover these 
parameters. The transfer function, H(Z) Is expressed as 


H(Z) 


2 _ 

D(Z) 


where fi(Z) Is the Infinite partial fraction expansion of 


D(Z) 


D(Zi 

F(Z) 


(2.26) 


(2.27) 


and the 4^ are the autoregressive coefficients. For practical applications 
this Infinite polynomial 6(Z) Is truncated after a finite number of terms. 

Graupe has derived error bands for this truncation. 

Since, In general, the user will have no aprlorl feel for idiere D(Z) should 
be truncated RAP solves for all coefficients for polynomials of order two 
thrott^ twenty. The solutions are obtained by recursive solution of the Yule- 
Walker equations which obtain the autoregressive coefficients as a linear func- 
tion of the autocorrelation coefficients as suggested by Box and Jenkins (Reference 

2 

2). The program also computes and plots the power transfer function |h(Z)| for 
the polynomial of ordar twenty. 
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3.0 TESTS FOR STATIOHARITY AND NORHALITY 

If an altimeter residual time series is stationary, this implies that the 
parameters of the underlying probability distribution are constants (l.e. not 
time varying) . RAF Includes the capability to assess stationarity and to test 
the hypothesis: "given that the series is stationary, then the underlying 

probability distribution is normal." 

Non stationarity would be indicated if the data contains significant 
secular (trend) or periodic conq>onents, or if the random component (noise) 
statistics change with time. To obtain a single assessment of stationarity 
RAP subdivides the time series into several segments and computes the mean and 
standard deviation for each segment. Overall "F" statistics will be computed 
to 1) assess the significance of differences in variability about the segment 
means, and 2) assess the significance of variability of the segment means 
about the overall mean. Because of the anticipated sample sizes, these statis- 
tics will be quite robust to departures of the underlying probability distribu- 
tion from normality. The power spectrum, described in Section 2.6 provides a 
more comprehensive assessment of stationarity. 

The data normality will be assessed by conq>uting chi square statistics 
which measure the difference between the sample distribution function and normal 
distribution function having the same mean and standard deviation. Kolmogorov- 
Smlmov two-sample test is used to assess the consistency of the underlying 
probability distribution function between segments. 


3.1 Chi Square Test for Distribution Normality 

To test for normality the sample histogram is tabulated and the following 
statistic conq>utcd (reference 3) 


2 

X 





(3.1) 


where r is the number of cells in the histogram, n^ is the number of sample 
values falling within the ith cell and e^ is theoretical expected number of 
sample values for the ith cell. The theoretical value, e^ is 
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evaluated aa 




n [ f(x)dx 


(3.2) 


where f(x) is the normal probability density function with mean and standard 
deviation estimated from the saiiq>le, is the lower cell boundary, is 

the upper cell boundary and n is the total sample size. The series approxi- 
mation used in SAP for i.he Integral of the normal distribution function was 
taken directly from reference 1. Since the mean and standard deviations are 

estimated from the sample the number of degrees of freedom associated with the 
2 

X statistic is r-2. Under the assugq>tlon that the underlying 8aiq>le dis- 

2 

trlbution is normal BAP computes the probability that x is less than or equal 
to the computed value. This is done by integrating the chi square density 
function as 


Prob (Xy _3 S 


I (frCOdt 

0 


(3.3) 


where ^(t) is the chi square density function with r-2 degrees of freedom. 
The method used in evaluating this Integral is described in reference 4. 


3.2 ”t” Statistic to Test for Zero Residual Mean 

If a saoq>le is drawn from a normal population the "t" statistic defined as 


^ X ** a 
sH n 


(3.4) 


follows a students "t" distribution with n-1 degrees of freedom where a is 
the hypothesized population mean. In RAP this statistic is computed for each 
column of data with a"0 to test the hypothesis that the column means are zero. 


3.3 P Test for Honoeeneity of Variances 

In order to test the hypothesis that a time series (l.e. a column of data) 



has constant variance the series Is partitioned into several segnents and the 
testing procedure reconnended in reference 2 is used. This procedure is 
suasiarlzed as follows: 

(1) Partition the time series into g segments and compute the trismied 
mean for each segment. The trlomted mean is obtained after deleting 
the largest ten percent and the smallest ten percent of the segpent 
data points. 

(2) Replace each data point by the absolute value of its deviation from 
the trimmed segment mean 


Hi " l*ij - *il 


(3.5) 


(3) C«q>ute the test statistic 


f n.(ii - z)^/(g-l) 
i-1 

f I^(*n ■ 

1-1 j-1 ^ 


(3.6) 


where 


'i • I *ij^"i 


i-1 j-l 


If the time series data consists of independent^ normally dietributed variates 
with constant variance then the test statistlCt W, will follow an "F" 
distribution with g-1, n-g degrees of freedom. Large values of W indicate 
that the variances are not constant. Results of monte carlo test’> summarized 
in reference 2 demonstrate that W is quite robust to departures from normalityi 
however, no analytic or numerical results were published to indicate what happens 
when the variates are autocorrelated. 


3.4 F Teat for Hoaogenelty of Mtans 


In order to test the hypothesis thet a time series (l.e. a column of data) 
has constant mean the series is partitioned into several segments and the 
following computational steps performed: 

(a) Confute the mean for each segiiient 
1 

X- ■ — y X 44 i«l,2***g where g - nund>er of segments (3.7) 

^ “i jtl 


(b) Compute the sum of squares about the mean for each segment 

^ v2 


asm - I (x - X ) 
^ j-1 ^ 


(3.8) 


(c) Confute the test statistic 


ilj^t^Xi-x)^/(g-l) 

\ ssm/(n-g) 
i -1 ^ 


(3.9) 


If the time series consists of independent normal variates with constant mean 
and variance then G follows an F distribution with g-1, n-g degrees of 
freedom. A large value of G and a "small" value of the W statistic (Eqn. 
( 3 . 6 )) would indicate the value of the mean is not constant. 


3.5 Kolmogorov-Smirnov Testing for Homogeneity of Sample Distribution Functions 

The homogeneity of the sample distribution function is tested, for each 
column of data, by partitioning the time series into several segments. Kolmogorov- 
Sedmov two-sample statistics are computed to compare the empirical distribution 
for each segment with the empirical distribution function obtained by grouping 
all of the remaining segments. Thus, if there are g segments there will be 
g test statistics computed. 

Following the procedure described in reference 5 the KolmogorovSmimov 
two'sample statistic is computed as follows: 
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Define ^2***^^ ^ sao^le of sice n and (y^^* y 2 ***ya^ ** 

the other sample of alae m. 

(a) The aaiqiles {Xj^, *2****n^ y2***^n^ *** sorted Into the 

ordered sets *(2)****(n)^ ^ ^^(1)* ^(2)‘“^(n)^ which 

are non decreasing sequences 

(b) The esq;>irical cumulative distribution functions and G^(y) 

are computed. For excon>le F (x) is defined by 

n 

0 for X < x^j 

F^(x) ■ V k/n for « x < h-l,2*«*n-3 (3.10) 

I 1 for X, V < X 
s, (o) 

(c) The maximum difference in absolute value between the two sample 
distribution functions is computed as 


D„ _ - max |F (x)-G (y)| 

Ulan n n * 

x,y 

The limiting distribution function L(z) for the test statistic 
given by (reference 5) 


(3.11) 



m.n 


is 


0 z ^ 0.27 

3 

{/1^/z) I expl-(2k-l)^ir^/8z^l + 
k-1 

Ej^(x); 0.27 < z < 1.0 

L(z) -( (3.12) 

r 2 2 

1-2 I (-1) exp (-2k V) + E»(z) 

k-1 


1.0 S z < 3.1 

^ 3* X ^ z 

vhtre: 

Ej^(x) t 6 (10“^^) when x < X 
E 2 (x) < 10“^® when x k 1 
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1.0 INTRODUCTION 

The SAMPFL Progran Is designed to filter atld subsai^le time series data for Input 
to the Residual Analysis Program (RAP). Three filtering options are Included. 

These are: 

• 6th order Buttervorth filter 


I 


• Adaptive polynomial filter 

• Fixed weight moving average filter 

The data may be low or high pass filtered. The three basic filters are low pass 
filters. The hljh pass output Is obtained by taking the difference (or residual) 
between the Input and the low pass filter output. 

Typically) a low pass filter would be used to remove high frequency noise 
prior to an analysis of the systematic components of the data. If a particularly 
long time series is to be analyzed to Heformine the low frequency spectral struc- 
ture of the data It is advisable to subsample the data OLfi@v low pass filtering 
to eliminate aliasing effects. 

A high pass filter would be used to remove the systematic or low frequency 
components of the time series prior to evaluation of the measurement noise 
statistics. 

SAMPFL contains logic for the extraction of data from an input tape between 
given start and stop tines. Outlying observations are deleted and replaced 
with a local average. The program also has the capability to either apply a 
fixed correction or entirely delete data between specified start and stop tlaws. 



2.0 FILTER DESCRIPTIONS 


Each of the three basic filters described in the following sections is a 
low pass filter, or smoother. The input is a discrete time sequence {y^) and 
the output is a discrete time sequence {x^}. Here the time sequence {x^} is 
smoother than {y^}. 

A discrete filter satisfies the following difference equation 

*1 ■ I »t n-t + ' \ *i-k 

If " 0 for all k then the filter is called a non recursive filter and 
output is simply a weighted average of the adjacent values of the input tiw 
series. If* in additon* the (g^) are constant* then the filter is called a fixed 
weight moving average filter. This type of filter is described in Section 2.1. 
If the {hy} are non sero the filter is called a recursive filter. Section 2.2 
describes a Butterworth filter for which the {g£,} and {hj^} are constant and 
Section 2.3 describes an adaptive polynomial filter which has time varying co> 


efficients. 


2.1 Plx»d W«lght Hovlog Averaga Filter 


:• 5 ' 


As defined in Section 2.0 {y^} Is the filter Input sequence end {x^} is 
the filter output sequence, With a fixed weight moving average filter each x^ 
is a linear combination of the 2n*fl surrounding values of the input sequence. 
Thus the filter output can be expressed as: 

-i “i 

tfhere {o^} is a fixed weight sequence. 

SAMFFLhas two moving average options. These are: 

(a) An arbitrary, externally determined set of weights {a^> stored In NXXX 
DATA are used; 

(b) The weights {o^} are computod so as to give a 2n+l point polynomial 
midpoint fit. 

For an (m-1)^^ order polynomial midpoint fit the regression equation is 
given by 

y, - + «i • 'o’ + •••■' Vi “i • V*'** ‘i 

where 1 ■ - n,..., 0...n and is assumed to be s zero mean random error. 
Since the data is assumed to be equlspaced in time, the time difference terms can 
be expressed as 


t. - t ■ iAt (?..4) 

1 o 

For equispacsd data the value of At is irrelevant in determining the weight 
sequence} therefore* in the following discussion it is assumed that At ■ 1.0. 

In matrix notation the least squares sstimate of the regression coefficients 
{6^ ... 6^,> is given by 

O Wrl, 

6 • (A^A)"^ A^ X (2*5) 

where v* - (y^ . . . y^ . . . y^ . ) contains 2wfl consecutive values of the ti'ilms 

^ t CtO 

series. The normal matrix A has 2n41 rows with the 1th row containing the 
partial derivatives of x^* as defined from eqvvition 2-3* with respect to the 


1*5 


regression coefficients. Thus the ith row of A is given by (for At ■ 1.) 


- [1 1 1^ ... l“*^l (2.6) 

For the filter aldpolnt 1 *• 0 and 

a^°^ - [1 0 0 ... 0] (2.7) 


The saoothed midpoint value is given by 

\ i - e. 


( 2 . 8 ) 


Substituting into (2.8) from (2.5), the expression for x^ in matrix nota- 
tion is given by 


- .«> (A^A)-l A^ i 


(2.9) 


Cooq>aring (2.9) with (2.2) which defines a moving average filter it is 
obvious that the weight sequence is given by 


- .«> (A^A)-l A^ 


( 2 . 10 ) 


Equation 2.10 is liq)lemented in SAMPFL to compute the weight sequence for 
an (m-l)th order polynomial midpoint fit spanning 2n't-l elements of the input 
tiM sequence. 

If a high pass filter is specified, the filter output is given by 

■ ’'i ■ 


(2.10a) 






2 . 2 Butterworth Filter 


Etpeating equation 2.1» a digital filter aatiaflea the difference equation 


-1 • J »i yi-* + j \ *i-k 


For a second order Buttervorth filter this speclallaea to 


*1 “ *1 ^1 ^^ 1-1 ^ 1 - 2 ^ *^1 * 1-1 * **2 * 1-2 


( 2 . 11 ) 


In SAMPFL a sixth order Buttervorth filter la Inplenented as three cascaded 
second order filters. The coefficients g^^* h^ and h2 are conputed for each 
stage as described In references 1 and 2 to give a ■axlnally flat frequency 
response in the paasband. 


The coefficients gj^, hj^ and h^ for the jth stage (J • 1»2,3) are given 


by 


0) 


g 


^ (I) ^ + a . (I) + 1 

a j * 


hj— 2gj ( 




( 2 . 12 ) 




kj - - gj ( 


u - a. u + 1 
_5 ) 

n / 




If u. la the desired cutoff frequency for the digital filter and At the 
d 


8aaq>llitg Interval then la given by 


u .At 

• ■ tan ( — g— ) (2.13) 

The coefficients (a^* j - 1»2»3) for a sixth order filter are given in 
reference 2 as 

Sj - .5176 

Sj - 1.4142 (2.14) 

Sj • 1.9318 
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2*3 Adaptive Polynomial Filter 


For a fixed length polynomial filter as described In Section 2*1« the length 
of the filter (2n+l) and the polynomial order (m-1) are arbitrarily specified. 

The adaptive polynomial filter atteiiq>ts to adjust the effective filter length to 
the data structure for an arbitrary polynomial order. The filter Is li^lemented 
as a Kalman filter i^lch adaptively adjusts a single state noise parameter. 

The following paragraphs describe the adaptive polynomial filter equations 
as Implemented In SAMPFL. In these discussions the subscripts of each of the 
vectors or matrices refer to first the data point at which the matrix or vector 
Is being evaluated, and second, the data point at which observation data was last 
Incorporated Into the estimate. Thus means the estimate prediction of 

£ at data point 1, based on observation data through data point 1-1. 
means the estimate of £ at data point 1-1, based on observation data through 
data point 1-1. evaluated at ■ S^.x/l-l* 

A particular component of a vector or matrix Is denoted by superscript. For 
example, first component of the vector l^/l-l* 

The state vector, ... contains regression coeffi- 

cients sladLlar to those defined In Section 2.1. 

The Incremental state model relating the filter state at tliK t^ to the 
state at t^^j^ “ ^1 defined by the equations 


( 0 ) 


( 0 ) 


* 1-1 


+ e 


( 1 ) 


1-1 


At + 


(m-1) 


* 1-1 


(At)-^ + q 


* q /Lt 

+ qAt 


e - 3 

^1 * 1-1 




(2.15) 


Here q is assumed to be a sero mean white noise process irlth unknown varlsnces 
Q. The measurement or Input, y^. Is related to the filter state es 


( 0 ) 


+ r^ 
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(2.16) 


where Is a zero mean white noise process with known variance 


The recursive filter update steps are as follows: 

(a) Update the filter state 

■^i/i-1 " *i.i-l -^-1/1-1 


(2.17) 


where the state transition matrix, is given by 


l.i-1 


1 At At 

01 0 

0 0 1 

• • • 

• • • 

0 0 0 


At 


m-1 


0 

0 


(b) Update the filter state covariance matrix 


*i/i-l “ *l,i-l ^i-l/i-1 *i,i-l * 


(2.18) 


(2.19) 


where Q is a scalar and r is a diagonal matrix having 

r. 


li 


(At)^“^ 1 - 1,2 ... m 


( 2 . 20 ) 


(c) Compute the residual variance 




(d) Conq>ute the filter gain matrix 


^ ^i 


i/i-1 


( 1 » 1 ) 1 


i/i-1 


( 1 , 2 ) 


L^i/1-1 


(l.m) 


( 2 . 21 ) 


( 2 . 22 ) 
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(f) Conpute the filter residual 


(g) Using the aeasureaent residual update the filter state estlaiate 

\/l ■ \/l-l ■'l .2.24) 

(h) Update the filter state covariance matrix to Include the measureoient 
information 

where is the m element row vector (1,0,0,...0) 

(i) Obtain an updated estimate of the scalar state noise parameter 


Vi " Vi-1 * 


(2.26) 


In SAHPFL, to Insure filter stability, is restricted to the range of 

values given by 


^in — ^i/1 - Snax 

where and are input parameters. 


(2.27) 
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